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ABSTRACT 

The long-term carbon cycle is vital for maintaining liquid water oceans on rocky 
planets due to the negative climate feedbacks involved in silicate weathering. Plate 
tectonics plays a crucial role in driving the long-term carbon cycle because it is respon¬ 
sible for CO 2 degassing at ridges and arcs, the return of CO 2 to the mantle through 
subduction, and supplying fresh, weatherable rock to the surface via uplift and orogeny. 
However, the presence of plate tectonics itself may depend on climate according to re¬ 
cent geodynamical studies showing that cool surface temperatures are important for 
maintaining vigorous plate tectonics. Using a simple carbon cycle model, I show that 
the negative climate feedbacks inherent in the long-term carbon cycle are uninhibited by 
climate’s effect on plate tectonics. Furthermore, initial atmospheric CO 2 conditions do 
not impact the final climate state reached when the carbon cycle comes to equilibrium, 
as long as liquid water is present and silicate weathering can occur. Thus an initially 
hot, CO 2 rich atmosphere does not prevent the development of a temperate climate 
and plate tectonics on a planet. However, globally supply-limited weathering does pre¬ 
vent the development of temperate climates on planets with small subaerial land areas 
and large total CO 2 budgets because supply-limited weathering lacks stabilizing climate 
feedbacks. Planets in the supply-limited regime may become inhospitable for life and 
could experience significant water loss. Supply-limited weathering is less likely on plate 
tectonic planets, because plate tectonics promotes high erosion rates and thus a greater 
supply of bedrock to the surface. 

Subject headings: astrobiology — planets and satellites; physical evolution — planets 
and satellites: terrestrial planets 


1. Introduction 

One of the biggest questions surrounding rocky extra-solar planets is whether they are poten¬ 
tially habitable for life. Water is thought to be necessary for the development of life, so a habitable 
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planet is typically defined as one that is able to sustain liquid water oceans (e.g. Kasting &: Catling 
2003). In order to be habitable a planet must lie within the habitable zone, the range of orbital 
distances where water can exist as a stable phase on a rocky planet’s surface (e.g. Hart 1978, 1979; 
Kasting et al. 1993; Franck et al. 2000), and it must have accreted enough water to produce oceans. 
However, lying within the habitable zone does not guarantee that a planet, even one that accreted 
a sufficient supply of H 2 O, will have liquid water oceans, or be able to maintain oceans for a sig¬ 
nificant portion of its history. The abundance of greenhouse gases in a planet’s atmosphere is also 
critical. A planet with a strong greenhouse effect could be hot enough to enter a moist greenhouse 
state, where dissociation of water in the upper atmosphere and hydrogen escape to space leads to 
rapid water loss (e.g. Kasting 1988; Abbot et al. 2012). Furthermore, stellar evolution exerts a 
major influence on planetary climate that must be balanced by atmospheric greenhouse gas abun¬ 
dances (Kasting 1989). Stars increase in luminosity as they age, so early in a planet’s history high 
greenhouse gas concentrations are needed to avoid a snowball climate (Sagan &: Mullen 1972; New¬ 
man & Rood 1977; Gough 1981), and late in a planet’s history low greenhouse gas concentrations 
are needed to keep surface temperatures below the moist greenhouse limit, where rapid water loss 
occurs (Kasting 1988). The same effects apply to planets at different orbital distances. Planets 
closer to the inner edge of the habitable zone need a lower greenhouse effect to maintain oceans, 
while those closer to the outer edge need a strong greenhouse effect to prevent global glaciation. 

A mechanism capable of regulating the strength of the atmospheric greenhouse effect, and 
thus counteracting the spatial and temporal variations in solar luminosity described above, exists 
via the long-term carbon cycle, which controls atmospheric CO 2 concentrations (e.g. Walker et al. 
1981; Kasting et al. 1993; Berner 2004). The long-term carbon cycle refers to the cycling of CO 2 
between the atmosphere and ocean, carbonate rocks on the seafloor, and the mantle. Weathering 
of silicate minerals on continents and in the oceanic crust draws CO 2 out of the atmosphere and 
ocean, depositing it on ocean plates in the form of carbonate rocks. CO 2 resides on the seafloor 
until it is subducted, where some portion of the carbon reaches the deep mantle, and the rest is 
returned to the atmosphere through metamorphic degassing and arc volcanism. The carbon in the 
mantle eventually degasses back to the atmosphere and ocean at mid-ocean ridges, completing the 
cycle. The long-term carbon cycle stabilizes planetary climate due to the temperature sensitivity 
of silicate weathering. At high temperatures weathering rates increase, and more rapid weathering 
draws CO 2 out of the atmosphere, cooling the climate. Meanwhile at low temperatures weathering 
rates decrease (or cease entirely in the case of an ice covered planet), and sluggish weathering allows 
volcanic degassing to build CO 2 back up in the atmosphere, warming the climate (e.g. Berner 2004). 
The same feedback buffers climate against changes in solar luminosity. When luminosity is low, 
slower weathering rates allow large quantities of CO 2 to buildup in the atmosphere, warming the 
climate, and when luminosity is high rapid weathering draws down atmospheric CO 2 , cooling the 
climate (e.g. Walker et al. 1981). As a result, an active global carbon cycle is thought to be essential 
for the long-term maintenance of liquid water oceans on rocky planets within the habitable zone 
(Kasting & Catling 2003). 
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Plate tectonics plays a vital role in the operation of the long-term carbon cycle. Plate tectonics 
drives volcanism at ridges and arcs, the major sources of atmospheric CO 2 , and facilitates silicate 
weathering, the primary sink of atmospheric CO 2 , by providing a continuous supply of fresh, 
weatherable rock (both on continents and on the seafloor) through orogeny and volcanic resurfacing. 
Therefore plate tectonics is often thought to be necessary for the operation of the long-term carbon 
cycle, and thus crucial for allowing planets throughout the habitable zone to sustain liquid water 
oceans over geologic timescales (Gonzalez et al. 2001; Kasting & Catling 2003) (though it is unknown 
whether some form of carbon cycling and climate stabilization can occur on non-plate-tectonic 
planets). Studies of carbon cycling and planetary habitability typically assume that plate tectonics 
and climate are independent; the presence of plate tectonics is imposed by the modeler, and the 
plate speed is often taken as a free parameter (Driscoll & Bercovici 2013), or as solely a function 
of mantle temperature (e.g. Tajika & Matsui 1992; Franck et al. 1999; Sleep & Zahnle 2001). 
However, recent geodynamical studies have found that climate and plate tectonics might be linked. 
Cool surface temperatures promote plate tectonics by suppressing grain-growth, allowing weak 
plate boundaries to form through grainsize reduction (Landuyt & Bercovici 2009; Foley et al. 2012; 
Bercovici & Ricard 2014), and by increasing mantle convective stresses, such that convective forces 
can more easily exceed the lithosphere’s intrinsic strength (Lenardic et al. 2008). Plate speed may 
even depend on surface temperature, as higher surface temperatures lead to stronger lithospheric 
shear zones via rapid grain-growth, and these stronger shear zones provide a larger resistance to 
plate motions (Foley &: Bercovici 2014; Foley et al. 2014). Although there is uncertainty over the 
applicability of these geodynamical studies to the tectonics of terrestrial planets, as the mechanism 
responsible for plate tectonics on Earth is still not fully known (Tackley 2000; Bercovici 2003), 
the implications of a fully coupled climate-mantle system (where the presence and vigor of plate 
tectonics depends on surface temperature) for planetary habitability are worth exploring. 

As stated in the previous paragraph, one of the key roles plate tectonics plays in the global 
carbon cycle is providing a supply of fresh rock at the surface. Without a sufficient supply of 
fresh rock continental weathering can enter a supply-limited regime. In the supply-limited regime 
the weathering reaction runs to completion in the regolith, and weathering can only continue by 
bringing unaltered bedrock into the weathering zone through physical erosion. As a result, supply- 
limited weathering depends on the physical erosion rate, rather than the kinetics of the weathering 
reaction (e.g. Stallard &: Edmond 1983; Edmond et al. 1995; Kump et al. 2000; Riebe et al. 2004; 
West et al. 2005; Hilley et al. 2010). If continental weathering on a planet becomes globally 
supply-limited, it would no longer depend on atmospheric CO 2 and surface temperature, and the 
climate stabilization provided by the long-term carbon cycle could be lost (West 2012). Given the 
important link between plate tectonics and the supply of weatherable rock, including a treatment 
for supply-limited weathering is crucial for studying the coupling between plate tectonics, climate, 
and planetary habitability. Furthermore, supply-limited weathering could be especially important 
on planets with small amounts of subaerial land (i.e. waterworlds), where the supply of fresh rock 
at the surface would naturally be low. Yet supply-limited weathering is often ignored in studies of 
exo-planet habitability, which, so far, have optimistically concluded that climate stabilization via 
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the long-term carbon cycle can still be maintained on planets where up to 99 % of the surface is 
covered by oceans (e.g. Abbot et al. 2012). 

Therefore in this study I test the ability of the long-term carbon cycle to regulate planetary 
climate, and thus allow hospitable surface conditions to persist over geologic time, for a fully 
coupled climate-mantle system that includes the supply limit to weathering. In particular, I assess 
whether the negative climate feedbacks inherent in the long-term carbon cycle are preserved when 
the plate speed depends on surface temperature as suggested by Foley & Bercovici (2014) and 
Foley et al. (2014), and the influence initial atmospheric CO 2 concentrations have on the final 
steady-state climate reached on a planet. Finally, I determine whether planets with small land 
areas can maintain a carbon cycle with negative climate feedbacks when supply-limited weathering 
is considered. This study only considers Earth-like planets, that is planets whose size and bulk 
composition are similar to the Earth’s, as these type of planets are most easily modeled based on our 
current understanding of climate, tectonics, and carbon cycling on Earth (the possible influences 
of planet size are discussed in §4.1.1). 

This study also focus solely on CO 2 cycling, which is thought to exert the leading order control 
on the climate evolution of rocky planets (e.g. Kasting & Catling 2003; Berner 2004). One reason 
for the prominence of CO 2 is that it is a major component of the gases released by terrestrial 
magmatism on Earth today, likely because the upper mantle became oxidized during, or soon after, 
core formation (e.g. Catling & Claire 2005; Wade & Wood 2005; Erost et al. 2008). Thus mantle 
degassing is a major source of atmospheric CO 2 , and without weathering to act as a sink for this 
CO 2 , Earth’s climate would be inhospitable to life. Carbon dioxide degassing has likely been active 
for the majority of Earths’ history, because core formation on the Earth was rapid, due to the high 
interior temperatures reached during accretion (Kleine et al. 2002). Geochemical data indicate 
that the upper mantle was oxidized by 3.9-3.5 Ga (Delano 2001; Frost & McCammon 2008), and 
possibly as early as ~ 4.4 Ga (Trail et al. 2011). Early oxidation of the mantle during core formation 
is thought to be a general feature of terrestrial planet formation, at least for planets as large or 
larger than the Earth (Halliday & Wood 2007), so CO 2 is likely to be a significant component of 
the atmospheres of rocky exo-planets as well. Although carbon dioxide is not the only greenhouse 
gas that can play a key role in planetary habitability, methane in particular may be important 
for solving the faint young sun problem (Pavlov et al. 2000; Haqq-Misra et al. 2008; Eeulner 2012; 
Charnay et al. 2013; Wolf &: Toon 2013), the ability of silicate weathering to regulate atmospheric 
CO 2 is likely to be crucial for the habitability of rocky planets. 

The paper is organized in the following manner: §2 describes the global carbon cycle model 
used in this study, with the weathering formulation that includes supply-limited weathering laid 
out in §2.2, and the dependence of plate speed on surface temperature described in §2.3; results 
are shown in §3, with the results of models testing the stability of the fully coupled climate-mantle 
system given in §3.1, and the influence of supply-limited weathering on climate stabilization for 
planets with small land areas shown in §3.2; finally, results are discussed in §4 and the main 
conclusions are summarized in §5. 
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2. Model Setup 


Global carbon cycle models of varying complexity have been used to study the long-term 
evolution of the Earth and other planets (e.g. Walker et al. 1981; Tajika & Matsui 1990, 1992; 
Franck et al. 1999; Sleep & Zahnle 2001; Sleep et al. 2001; Abbot et al. 2012; Driscoll & Bercovici 
2013), as well as for more detailed studies of Neoproterozoic and Phanerozoic climate evolution (e.g. 
Berner et al. 1983; Volk 1987; Berner 1994, 2004; Tajika 1998; Mills et al. 2011). In this study a 
simple model is used to focus on the first order effects of both climate-tectonic coupling and exposed 
land area on climate stabilization via the global carbon cycle. Carbon is assumed to cycle between 
four reservoirs: the ocean crust (or plate) (Rp), the mantle (Rman), the atmosphere (Ratm), and the 
ocean (Roc) (Figure 1). A significant fraction of oceanic carbonates reside on continental shelves at 
the present day (e.g. Ronov &: Yaroshevsky 1969), and thus a separate continental reservoir is often 
included in global carbon cycle models (e.g. Tajika &: Matsui 1990; Sleep & Zahnle 2001). Here all 
carbonates are assumed to form on the seafloor, i.e. they are part of the ocean crust reservoir, as 
in Driscoll &: Bercovici (2013). The separate continental reservoir is neglected because much of this 
study focuses on planets with small land fractions, where the large majority of carbonate formation 
would occur on the seafloor. Furthermore, neglecting the continental reservoir does not significantly 
alter the first order climate feedbacks or climate-tectonic coupling, because the functional form of 
the degassing function for continental carbon is similar to the mantle degassing function, in that 
both depend linearly on plate speed (Sleep & Zahnle 2001). 

Carbon cycles between the four reservoirs in this model as follows (Figure 1): Exposed rock at 
the surface reacts with atmospheric CO 2 to form bicarbonate, magnesium, and calcium ions that 
travel to the oceans via rivers and groundwater. Once in the ocean, these ions form carbonate 
minerals on the seafloor, both in the form of sediments and via hydrothermal alteration of basalt 
(CO 2 dissolved in the oceans can also react directly with basalt and form carbonate minerals, 
regardless of whether silicate weathering is active) (e.g. Berner 2004). Thus weathering acts as 
a net flux of carbon from the atmosphere and ocean to the plate reservoir. Carbon leaves the 
plate reservoir when it is subducted at trenches. Here, a fraction of the carbonate minerals will 
devolatilize, returning carbon to the atmosphere via arc volcanism, and the remainder will be 
subducted into the deep mantle. Once in the mantle, carbon is mixed throughout by convection, 
and degasses to the atmosphere and ocean at mid-ocean ridges. Return of carbon to the atmosphere 
and ocean via mantle degassing closes the long-term carbon cycle. 


The carbon cycle is modeled using the following mass balance equations (e.g. Tajika &: Matsui 
1992; Sleep & Zahnle 2001; Driscoll Sz Bercovici 2013): 
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where F^eather (the continental or terrestrial weathering flux) is the flux of carbon from the atmo¬ 
sphere to the plate resulting from silicate weathering on land, Fgfw (the seafloor weathering flux) 
is the flux of carbon from the ocean to the plate resulting from hydrothermal alteration of basalt 
by CO 2 dissolved in seawater, Fg^h (the subduction flux) is the flux of carbon from the plate into 
the mantle, F^egas (the degassing flux) is the flux of carbon from the mantle to the atmosphere and 
oceans via mid-ocean ridge volcanism, and Fare (the arc volcanic flux; Fare = fFgubj where / is the 
fraction of subducted carbon that is degassed at arcs) is the flux of carbon to the atmosphere from 
devolatilization of the subducting slab (see Table 1 for a list of all model variables). The continen¬ 
tal weathering flux is divided by two because half of the carbon initially drawdown by weathering 
on land is rereleased to the atmosphere when carbonates form on the seafloor (e.g. Berner et al. 
1983). In other words, Fu,eather/‘^ represents the net flux of CO 2 from the atmosphere and ocean 
reservoirs to the plate reservoir via continental weathering. The seafloor weathering flux only rep¬ 
resents the alteration of basalt by CO 2 dissolved in the oceans; formation of carbonates in basalt 
by calcium and bicarbonate ions derived from silicate weathering on land are included in Fa,eather- 
Thus, equation (1) describes how carbon is added to the plate reservoir through continental and 
seafloor weathering, and leaves via subduction, equation (2) tracks the addition of carbon to the 
mantle reservoir via deep subduction of carbonates (i.e. (1 — f)Fsub), and its removal by ridge 
degassing, and equation (3) describes how arc volcanism and mantle degassing add carbon to the 
atmosphere and ocean, while silicate weathering removes it. This model also assumes that carbon 
is instantaneously mixed throughout the mantle after being subducted. This is a simplification 
because mantle mixing is not instantaneous, and can take up to 1 Gyr for whole mantle mixing in 
the present day Earth (e.g. Olson et al. 1984; Hoffman &: McKenzie 1985; Christensen 1989; Kellogg 
&: Turcotte 1990). However, as mantle mixing with realistic geometries and rheological effects is 
still a major topic of research (Tackley 2007), incorporating the influence of a finite mixing time 
for mantle carbon is left for future studies. 

As in Sleep &: Zahnle (2001), the atmosphere and ocean reservoirs are grouped together. The 
partitioning of carbon between the atmosphere and ocean is controlled by the solubility of CO 2 in 
the oceans. These two reservoirs equilibrate on a very short timescale compared to the timescale 
over which the other reservoirs evolve. Thus, equilibration between the atmosphere and ocean is 
assumed to take place instantaneously, and the combined atmosphere plus ocean carbon reservoir 
is partitioned into Ratm and Rqc using Henry’s Law at each timestep. Henry’s Law is given by 
Pc 02 — keXc, where Pco 2 is the partial pressure of atmospheric CO 2 , ke is the solubility and Xe 
is the mole fraction of CO 2 in the oceans; Xe = Roc/{MH2O + Roc) where MH2O is the number 
of moles of H 2 O in an ocean mass of water. The temperature dependence of the solubility, kc, is 
ignored in this study for simplicity; a kc representative of warm greenhouse conditions, where less 
CO 2 can be dissolved in the oceans, is chosen and used in the majority of the models (see Table 
2 for a list of all model parameters and their assumed baseline value, and Table 3 for a list of all 
physical constants used in this study). Variations in kc are considered in §3.2.1 and found to have 
no significant influence on the results. The atmospheric CO 2 reservoir can be related to the partial 


pressure of atmospheric CO 2 {PCO 2 ) through 
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Pc02 


Patm ^C 02 9 


■dearth 


( 4 ) 


where rhc 02 molar mass of CO 2 , g is acceleration due to gravity, and A^arth is the surface 

area of the Earth. This relation is only an approximation, because relating the moles of atmo¬ 
spheric CO 2 to Pco 2 depends on the composition of the rest of the atmosphere. If CO 2 is the 
dominant component, then equation (4) is exact; if N 2 or H 2 O are the dominant components in the 
atmosphere, then equation (4) will overestimate Pco 2 by approximately a factor of 2. However, 
the impact of this error on the overall results of the study is negligible; test cases assuming that 
CO 2 is a minor atmospheric component were found to be nearly identical to the results presented 
in the remainder of this paper, all of which use equation (4). 


2.1. Subduction, Arc, Degassing, and Seafloor Weathering Fluxes 


The subduction flux, Fgub, is given by the product of the area density of carbon on the plate, 
Rp/Ap (where Ap is the area of oceanic plates), the plate speed, v, and the length of trenches, L, 


RpvL 

Ap 


( 5 ) 


The arc flux is then Fare = fPsub: and the flux of carbon into the deep mantle is (1 — f)Fsub- The 
fraction of subducting carbon that degasses through arc volcanism on the present day Earth is not 
well constrained, with estimates typically ranging from 25 % to 70 % (e.g. Sleep & Zahnle 2001; 
Dasgupta & Hirschmann 2010; Ague &: Nicolescu 2014), though a vary small / < 0.1 can not be 
ruled out (Kelemen &: Manning 2015). Thus / = 0.5 is chosen as a baseline; variations in / are 
considered in §3.2.1, and do not signihcantly influence the main results of this study (see also §4.2). 


The degassing flux is given by the flux of upwelling mantle into the melting region beneath 
the ridge, multiplied by the concentration of carbon in the mantle and the fraction of this carbon 
that degasses (Tajika &: Matsui 1992). The flux of mantle into the melting zone must balance the 
flux of mantle leaving this region, and is thus given by 2vLdmeit, where dmeit is the depth where 
mid-ocean ridge melting begins and the length of ridges is assumed to be equal to the length of 
trenches. The degassing flux is then 


Pdegas — fd 


R„ 


14 


2vLdraelt 


( 6 ) 


where fd is the fraction of upwelling mantle that degasses and Vmantie is the volume of the mantle. 
The current day length of mid-ocean ridges is ~ 6 x lO'^ km (e.g. Fowler 2005); this value is assumed 
to be constant and used for all models. The degassing fraction, fd, is fixed to reproduce the present 
day Earth’s degassing flux and atmospheric CO 2 content, which gives fd ~ 0.32, similar to the 
estimate of Tajika & Matsui (1992). 
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Fig. 1.— Schematic diagram of the global carbon cycle model used in this study. 


Table 1: Table of variables. 


Symbol 

Definition and units 

Equation 

Rp 

Ocean plate carbon reservoir (mol) 

(1) 

Rman 

Mantle carbon reservoir (mol) 

(2) 

Ratm 

Atmosphere carbon reservoir (mol) 

(3) 

Roc 

Ocean carbon reservoir (mol) 

(3) 

Rsub 

Subduction flux (mol Ma“^) 

(5) 

Rdegas 

Degassing flux (mol Ma“^) 

(6) 

Rweather 

Continental weathering flux (mol Ma“^) 

(14) 

7? 

^ Ws 

Supply-limited weathering flux (mol Ma“^) 

(10) 

7? 

Kinetically-limited weathering flux (mol Ma“^) 

(11) 

Rsfw 

Seafloor weathering flux (mol Ma“^) 

(7) 

PC02 

Partial pressure of atmospheric CO 2 (Pa) 

(4) 

Ap 

Area of ocean plates {AEarthfiand, km^) 

(5) 

hand 

Land fraction 

(9) 

Psat 

Saturation vapor pressure (Pa) 

(12) 

T 

Surface temperature (K) 

(15) 

Te 

Effective temperature (K) 

(16) 

s 

Solar irradiance (W m“^) 

(16) 

V 

Plate speed (cm yr“^ or m Ma“^) 

(17) 

Ctot 

Total planetary CO 2 budget (mol) 

(18) 
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The seafloor weathering flux, in particular whether it has a strong climate feedback or not, 
is poorly understood. Drill cores of oceanic crust show that low temperature (0-60° C) off-axis 
hydrothermal alteration of basalt is a significant CO 2 sink (Staudigel et al. 1989; Alt & Teagle 
1999; Gillis & Coogan 2011) and laboratory measurements indicate that the rate of CO 2 uptake 
via basalt alteration increases with both increasing temperature and Pco 2 (Brady &: Gislason 1997). 
However, it is not clear if the temperature of the seawater circulating through crustal basalt during 
hydrothermal alteration is dictated by the geothermal heat flow (Alt & Teagle 1999) or by the 
ocean bottom water temperature (Gillis & Googan 2011; Googan & Gillis 2013; Coogan &: Dosso 
2015), which is related to climate. Furthermore, seafloor weathering can also become “supply- 
limited,” if all of the basalt accessible to hydrothermal fluids is completely carbonated (Sleep et al. 
2001); in this case only the creation of new seafloor at mid-ocean ridges allows further seafloor 
weathering. Thus even with a direct temperature feedback there are limits to the amount of CO 2 
that can be sequestered in the ocean crust. In this study both the direct temperature dependence 
of seafloor weathering and its supply-limit are neglected due to the large uncertainties associated 
with these effects; this also results in a simpler model that can be more completely understood 
and analyzed. In §3.2.3 I show that including a direct temperature feedback and a supply-limit to 
seafloor weathering does not significantly change the results of this study, therefore justifying the 
exclusion of these two effects, as long as complete basalt carbonation cannot penetrate to great 
depths within the crust. Physically, the simplified seafloor weathering flux used in this study can 
be interpreted as assuming that crustal heat flow sets the temperature of the water reacting with 
basaltic crust (Abbot et al. 2012). The seafloor weathering flux therefore follows after Sleep & 
Zahnle (2001) and Mills et al. (2014) as. 


Fsfw — Psfw 




(7) 


where stars denote present day quantities, a describes the dependence of basalt carbonation on 
atmospheric GO 2 , and {vjv*) represents the effect of spreading rate, with a constant ridge length, L. 
Spreading rate sets the rate at which fresh, weatherable basalt is created, and the Pco 2 dependence 
is based on Brady & Gislason (1997), who find a « 0.23; this is rounded up to 0.25 for the baseline 
value of a used in this study. The current day seafloor weathering flux, is set to 1.75 x 10^^ 

mol/yr (Mills et al. 2014). 


2.2. Continental Weathering Flnx 

The continental weathering flux typically used in global carbon cycle models is based on the 
kinetics of the weathering reaction between GO 2 and silicate rocks, with parameterizations for 
additional effects, such as increases in weathering rates brought about by increases in runoff or 
tectonic activity (e.g. Walker et al. 1981; Tajika &: Matsui 1992; Sleep & Zahnle 2001; Berner 2004; 
Driscoll & Bercovici 2013). However, the weathering flux will only be determined by the reaction 
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kinetics when the weathering reaction is the rate limiting step (i.e. when weathering is kinetically- 
limited). When weathering rates are high, or physical erosion rates are low, weathering can be 
limited by the supply of fresh rock to the surface (e.g. Stallard &: Edmond 1983; Edmond et ah 
1995; Kump et al. 2000; Riebe et al. 2004; West et al. 2005). A more general weathering function, 
that incorporates supply-limited weathering, is given by (Gabet &: Mudd 2009; Hilley et al. 2010; 
West 2012), 


W = excc 


1 — exp ( ^ j 


( 8 ) 


where W is weathering rate (with dimensions of mass area“^ time“^), e is the total denudation 
rate (mass area“^ time“^), Xcc is the fraction of reactable cations in the continental crust and is 
unitless, R is the rate of the silicate weathering reaction (time“^), and z is the effective depth 
of the weathering zone (mass area“^). The weathering rate, W, can be converted into the global 
weathering flux, Fweather, by integrating over the area of exposed land, and dividing by the molar 
mass of reactable cations in the bedrock. Thus, 


Fweather — 


-dearth fland 


rrir 


(9) 


where fiand is the fraction of exposed land (i.e. the area of exposed land divided by the surface area 
of the Earth), and rhcc is the average molar mass of reactable elements in the continental crust. 
Combining equations (8) and (9), 


(1 _ ,xp (-H?)) = (l - «p (10) 

where E is the physical erosion rate (length time“^), pr is the density of the regolith, and = 
{AearthfiandExccPr)/Flee is the Supply limit to weathering. Physical erosion rates vary by a wide 
range on Earth, depending on relief, climate, lithology, and other factors. In this study, a full model 
of how erosion varies in response to these factors is not attempted. Instead, a maximum erosion rate, 
Emax, representing the upper bound on the globally averaged erosion rate on a planet, is chosen. 
This maximum erosion rate, along with the area of exposed land, then sets the global supply limit 
to weathering, E^^. Emax is poorly constrained, so a wide range of maximum erosion rates are 
tested in §3.2.1, and the likely dependence of Emax on tectonic activity, and the implications this 
has for weathering and climate regulation on stagnant lid planets, is discussed in §4.4. 


The reaction rate, R, is determined from the typical weathering flux equations used in previous 
global carbon cycle models (e.g. Walker et al. 1981; Tajika & Matsui 1992; Sleep & Zahnle 2001; 
Berner 2004; Driscoll & Bercovici 2013); in this way the total weathering flux will mimic previous 
models when weathering is kinetically-limited, and then plateau when weathering rates are high 
enough, or erosion rates low enough, to reach the supply limit. Following Driscoll &: Bercovici 
(2013), the weathering flux for kinetically-limited weathering is 
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where is the present day weathering flux [F^ ss 12 x 10^^ mol/yr (Gaillardet et al. 1999)), Psat is 
the saturation vapor pressure, /3 and a are constants, Ea is the activation energy of the weathering 
reaction, Rg is the universal gas constant, T is the surface temperature, and starred quantities 
represent present day values. The weathering flux increases with increasing Pco 2 because 

both factors increase the rate of the weathering reaction (Berner 2004), and the saturation vapor 
pressure term models the variation in runoff as surface temperature changes (Driscoll &: Bercovici 
2013). The saturation vapor pressure is (e.g. Kasting et al. 1984; Nakajima et al. 1992) 


Psat — PsatO exp 


'^wFw ( 1 

~^\f 


1 

FsatO 


( 12 ) 


where Pgato is the reference saturation vapor pressure, at reference temperature Tgato, mw is the 
molar mass of water, and is the latent heat of water (see Tables 2 &: 3). 

Assuming that the kinetically-limited global weathering flux is equal to the reaction rate, R, 
multiplied by the total amount of reactable cations within the weatherable bedrock, 


F — R 


hw dearth flandPrXcc 


rUr. 


= R 






rrir 


(13) 


where is the thickness of the weathering zone, which can be replaced by zjpr (West 2012). 
Combining equations (10), (11), and (13), the definition of and using e = EmaxPr, 


F 


weather 


= F, 


1 — exp 


Fwfland f PC02 \ 
P^sfland \^C02 / 




(14) 


Inverting for the reaction rate using a weathering function normalized to the present day Earth 
(equation (13)) implicitly assumes that weathering on the present day Earth is entirely kinetically 
controlled. While the breakdown of the present day global weathering flux between supply-limited 
and kinetically-limited weathering is not well constrained, many estimates find that kinetically- 
limited weathering is the dominant contributor (Kump et al. 2000; West et al. 2005). Thus as¬ 
suming Earth’s present day weathering flux is kinetically controlled is a reasonable first order 
approximation. 

The weathering flux, F^f,ather, is shown in Figure 2, plotted against F^^., i.e. the typical weath¬ 
ering flux used in global carbon cycle models that do not treat supply-limited weathering. F^^ather 
is identical to at low Pco 2 (corresponding to low surface temperatures and low weathering 
rates), and diverges from F^^. at high Pco 2 when weathering becomes supply-limited; here F^^ather 
plateaus at F^^, the global supply limit to weathering, determined by the maximum erosion rate 
and area of exposed land at the surface. 
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Atmospheric CO2 [Pa] 

Fig. 2.— Comparison between the total weathering flux, including supply-limited weathering, 
Fweather (solid line), and the kinetically-limited weathering flux, (dashed line). Surface tem¬ 
perature depends on Pco 2 based on equation (15) and saturation vapor pressure depends on surface 
temperature as given by equation (12). A maximum erosion rate of 10 mm/yr and Earth-like land 
fraction of 0.3 were assumed, and the other parameters are listed in Tables 2 &: 3. 









- 13 - 


2.3. Climate and Plate Tectonic Models 


To relate atmospheric CO 2 to surface temperature, a simple parameterization from Walker 
et al. (1981) is used: 


T = T* + 2{Te - T*) + 4.6 


Pc02 
p, 


0.346 


-4.6 


(15) 


C02 


where T* = 285 K is the present day surface temperature, Tg is the effective temperature, and 
T* = 254 K is the present day effective temperature (Figure 3a). Note that this parameterization 
also includes the contribution of water vapor to greenhouse warming, assuming that H 2 O is always 
saturated in the atmosphere. The effective temperature is related to the absorbed solar radiation 
as 


T, = 


5(1 


Aa 


(16) 


where 5 is the solar irradiance, A is the albedo, and a is the Stefan-Boltzmann constant. In 
this study albedo will be held constant for simplicity, so ice-albedo feedbacks, or the currently 
poorly understood feedbacks involving clouds (e.g. Leconte et al. 2013; Wolf &: Toon 2014) are not 
included. This parameterization provides a good first order approximation to the results of more 
sophisticated radiative-convective models for an Earth-like planet (e.g. Kasting & Ackerman 1986). 
As discussed in §4.2, using a more advanced climate model may change some of the details of the 
results, particularly the temperatures calculated at high CO 2 levels where equation (15) deviates the 
most from the models of Kasting & Ackerman (1986), but would not change the overall conclusions 
of this study. 


Although temperatures reach ~ 600 K at Pco 2 ~ liquid water is still stable at the 

surface. A true runaway greenhouse, where liquid water cannot exist as a stable phase on the 
surface, is thought to only be possible through increases in insolation, rather than through the 
increase in atmospheric opacity brought about by elevated CO 2 levels (Nakajima et al. 1992). 
However, above ~ 350 K the atmosphere would be in a moist greenhouse state, where, although 
liquid water is stable at the surface, high mixing ratios of H 2 O in the stratosphere lead to rapid 
water loss to space (Abbot et al. 2012). Such water loss is not modeled in this study; rather I look 
at the expected climate and tectonic states for terrestrial planets with varying land areas, CO 2 
budgets, and incoming solar fluxes. However, the implications of the modeled climate states for 
water loss and habitability are discussed in §3.2.2 and §4. 


As discussed in §1, surface temperature can exert a major influence over whether plate tectonics 
takes place on a planet (Lenardic et al. 2008; Landuyt & Bercovici 2009; Foley et al. 2012), and 
influences the speed of plate tectonics when it is present (Foley & Bercovici 2014; Foley et al. 2014). 
To include the effect of climate on plate tectonics, plate speed is assumed to depend on surface 
temperature using a linear approximation to the scaling laws of Foley &: Bercovici (2014) 


V = 8.059 - O.OllT, 


(17) 




Surface Temperature [K] 
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^ C02 [1^^] 



Surface Temperature [K] 


Fig. 3.— Surface temperature as a function of Pco 2 for a present day solar luminosity and albedo 
(a), and plate speed as a function of surface temperature (b). Both the plate speed parameterization 
used in this study, equation (17), and the plate speed-surface temperature curve calculated from 
the scaling laws of Foley &: Bercovici (2014) are shown in (b). 
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where high surface temperatures lead to slower plates due to increased grain-growth rates in the 
lithosphere, which cause effectively stronger plate boundaries. Equation (17) gives v in units of 
cm yr“^. Mantle temperature also exerts an influence on plate speed. However, the influence of 
mantle temperature on plate speed is neglected in this study because it is weaker than the influence 
of surface temperature (Foley et al. 2014). Furthermore, ignoring the weaker influence of mantle 
temperature on plate speed also allows this study to focus solely on the coupling between plate 
tectonics, carbon cycling, and climate. Figure 3b shows that the plate speed parameterization, 
equation (17), is an accurate approximation of the full scaling laws from Foley & Bercovici (2014) 
up to Ri 600 K, corresponding to the highest surface temperatures explored here. 


3. Results 

In this section, I first explore the influence the full coupling between plate tectonics and climate 
has on the development of habitable conditions (§3.1), by: testing whether the ability of the long¬ 
term carbon cycle to buffer climate against changes in solar insolation is maintained when plate 
speed depends on surface temperature (§3.1.1); and by assessing whether habitable climates can 
develop from arbitrary initial surface temperature and atmospheric CO 2 conditions (§3.1.2), or 
whether different initial conditions prevent temperate climates from forming. Next, I examine the 
role land area and planetary CO 2 budget play in the development and maintenance of habitable 
climates, and map out the conditions necessary for stable, temperate climates to exist on rocky 
planets (§3.2). Throughout this section, all parameters are held constant (see Tables 2 & 3), save 
for plate speed which is allowed to vary with surface temperature, so that the coupling between 
climate and plate tectonics, and the influence of land area and total CO 2 budget, can each be 
studied in isolation. The effect of varying key model parameters is explored in §3.2.1, and the 
influence of including both the temperature dependence of seafloor weathering and its supply-limit 
is tested in §3.2.3. 


3.1. Carbon Cycling and Climate Feedbacks with a Surface 
Temperature-Dependent Plate Speed 

3.1.1. Stability of Climate in Response to Changes in Insolation 

As outlined in §1, one of the most important roles the long-term carbon cycle plays in planetary 
habitability is stabilizing climate against changes in solar luminosity. In this section, the efficacy 
of the carbon cycle climate buffer, when the influence of surface temperature on plate speed is 
included, is assessed. The steady-state solution to the global carbon cycle model (equations (1)- 
(3)) as S is varied is shown in Figure 4. Solar luminosities ranging from 0.8S*, appropriate for the 
early Earth (Gough 1981), to 1.25*, where S* = 1360 Wm“^ is the solar constant, are considered. 
Solving the full model in steady-state assumes that the global carbon cycle reaches steady-state 
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Table 2: 

Table of parameters. 



Symbol 

Definition 

Baseline Value 

Equation 

kc 

Solubility of CO 2 in seawater 

10^ Pa (derived) 

above (4) 

Mh20 

Moles of H 2 O in one ocean mass 

7.6 X 10^^ mol (derived) 

above (4) 

L 

Length of trenches 

6 X 10*^ km (F05) 

(5) 

f 

Fraction of subducted carbon that degasses 

0.5 (A14) 

(5) 

fd 

Fraction of upwelling mantle that degasses 

0.32 (T92) 

(6) 

dmelt 

Depth of melting beneath ridges 

70 km (K06) 

(6) 

^ sfw 

Present day seafloor weathering flux 

1.75 X 10^® mol Ma“^ (M14) 

(7) 

a 

Pco 2 exponent for seafloor weathering 

0.25 (B97) 

(7) 

* 

V 

Present day plate speed 

5 cm yr“^ (K06) 

(7) 

Xcc 

Fraction of Mg, Ca, K, and Na in continental crust 

0.08 (W12) 

(8) 

fhcc 

Average molar mass of Mg, Ca, K, and Na 

32 g mol“^ (derived) 

(9) 

E 

^max 

Maximum erosion rate 

10 mm yr“^ (W13) 

(10) 

Pr 

Regolith density 

2500 kg m-3 (W12) 

(10) 

P 

Pco 2 exponent for silicate weathering 

0.55 (D13) 

(11) 

a 

Psat exponent for silicate weathering 

0.3 (D13) 

(11) 

Ea 

Activation energy for silicate weathering 

42 kJ mol~^ (B91) 

(11) 

/?* 

W 

Present day weathering flux 

12 X 10^® mol Ma"^ (G99) 

(11) 

P* 

^C02 

Present day atmospheric CO 2 

33 Pa (K84) 

(11) 


Present day surface temperature 

285 K (W81) 

(11) 

land 

Present day land fraction 

0.3 (A12) 

(11) 

PsatO 

Reference saturation vapor pressure 

610 Pa (K84) 

(12) 

Lw 

Latent heat of water 

2469 J g-^ (K84) 

(12) 

TsatO 

Reference temperature 

273 K (K84) 

(12) 

e 

Present day effective temperature 

254 K (W81) 

(15) 

A 

Albedo 

0.31 (derived) 

(16) 


Key for citations: A12=(Abbot et al. 2012), A14=(Ague &: Nicolescu 2014), B91=(Brady 1991), 
B97=(Brady & Gislason 1997), D13=(Driscoll & Bercovici 2013), F05=(Fowler 2005), G99=(Gaillardet 
et al. 1999), K84=(Kasting et al. 1984), K06=(Korenaga 2006), M14=(Mills et al. 2014), T92=(Tajika 
&: Matsui 1992), TS=(Turcotte & Schubert 2002), W81=(Walker et al. 1981), W12=(West 2012), 
W13=(Willenbring et al. 2013). 





- 17 - 


on a timescale shorter than changes in solar luminosity; in reality the two timescales are similar, 
both on the order of ~ 1 Gyr. A different approach would be to assume that the plate and 
mantle reservoirs are fixed, and solve for the atmospheric reservoir in steady-state (i.e. setting 
Fsfw + Fweather/‘^ = f Fsub + Fdegas)- Both approaches were found to give nearly identical results. 

Similar to other studies on climate stabilization and the long-term carbon cycle (e.g. Walker 
et al. 1981), the temperature dependence of kinetically controlled weathering allows higher atmo¬ 
spheric CO 2 levels to build up at low solar luminosity, and lowers atmospheric CO 2 levels through 
higher weathering rates when solar luminosity increases (Figure 4B). As a result, the changes in 
surface temperature brought about by variations in solar luminosity are less pronounced than the 
case where atmospheric CO 2 is fixed (Figure 4A). Also consistent with previous studies (e.g. Sleep 
&: Zahnle 2001), the carbon cycle is not able to produce atmospheric CO 2 levels high enough to 
keep surface temperatures above 273 K at solar luminosities lower than ~ 0.95*. Other greenhouse 
gases, such as methane, may be necessary to prevent global glaciation. However, recent results 
from three-dimensional global climate models have shown that low latitudes can remain ice free 
at globally averaged surface temperatures of 260 K (Charnay et al. 2013; Wolf & Toon 2013), and 
possibly even lower (Abbot et al. 2011), so the ability of the global carbon cycle to maintain liquid 
water oceans at low solar luminosities may be stronger than previously thought. 

Figure 4 shows that the dependence of plate speed on surface temperature has almost no 
influence on the effectiveness of carbon cycle induced climate stabilization. The weathering feedback 
is by far the dominant factor, and thus the fully coupled climate-plate tectonic system is stable 
to changes in solar luminosity. In fact, the influence of surface temperature on plate speed has 
a very small additional stabilizing effect on long-term climate. Plate speed increases at cooler 
temperatures, so lower insolation causes degassing to increase (both at volcanic arcs and mid-ocean 
ridges), in addition to lowering the weathering rate. A complimentary effect occurs at higher 
insolation, where the warmer surface temperature slows plates and decreases degassing. However, 
as plate speed only changes by ss 0.2 cm/yr per 20 degree change in surface temperature, while the 
weathering flux changes by an order of magnitude, the feedbacks inherent in silicate weathering are 
by far the dominant factor in long-term climate stabilization. 


3.1.2. Influence of Initial Conditions 

The results in §3.1.1 show that habitable climates can be maintained even when plate speed 
depends on surface temperature. However, they do not demonstrate whether habitable climates 
develop in the first place. Fully time-dependent evolution models are needed to determine whether 
different initial conditions, especially different initial atmospheric CO 2 concentrations, influence 
the final surface temperatures reached when the global carbon cycle comes to steady-state. In 
this section, the full time-dependent form of the carbon cycle model, equations (l)-(3), is solved 
using three different initial conditions: an initially hot case, where the entire CO 2 budget of the 
planet is initially in the atmosphere and ocean, an initially cold case where the entire CO 2 budget 
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Table 3: Table of physical constants. 


Symbol 

Definition 

Value 

Equation 

rnc02 

Molar mass of CO 2 

44 g mol~^ (derived) 

(4) 

9 

Acceleration due to gravity 

9.8 m s-2 (TS) 

(4) 

^mantle 

Volume of the mantle 

9.1 X km^ (TS) 

(6) 

-^Earth 

Surface area of Earth 

5.1 X 10^^ km^ (TS) 

(4) 

Rg 

Gas constant 

8.314 J K-i rnol-i (N92) 

(11) 

niw 

Molar mass of water 

18 g mol~^ (derived) 

(12) 

a 

Stefan-Boltzmann constant 

5.67 X 10-® W m-2 K-^ (N92) 

(16) 

s* 

Solar constant 

1360 W m“^ (dPlO) 

§3.1.1 


Key for citations: dP10=(de Pater & Lissauer 2010), N92=(Nakajima et al. 1992), 
TS=(Turcotte Sz Schubert 2002). 




Solar Luminosity (S/S*) Solar Luminosity (S/S*) 


Fig. 4.— Surface temperature (a) and atmospheric CO 2 (b) as functions of insolation normalized by 
the solar constant. Surface temperature and atmospheric CO 2 are determined by solving equations 
(l)-(3) in steady-state for each S/S*. The black line assumes CO 2 is constant and equal to the 
present day value {Pco 2 = -^ 002 )’ shows the results from the carbon cycle model with 

a plate speed that depends on surface temperature, and the red line shows the results when plate 
speed is hxed at 5 cm/yr. Note that the blue and red lines lie almost on top of each other, and are 
therefore difficult to distinguish. 
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initially resides in the mantle, and an intermediary case where half the CO 2 starts in the atmosphere 
and ocean, and half in the mantle. An Earth-like planet is assumed, with a total CO 2 budget of 
Ctot = 2.5 X 10^^ moles (Sleep & Zahnle 2001), and a constant land fraction of 0.3. Both the 
thermal evolution of the mantle and the evolution of solar luminosity are ignored in this section so 
that the dynamics of the coupled climate-plate tectonic system can be studied in isolation. 

Figures 5, 6, and 7 demonstrate that all three initial conditions reach the same steady state 
after ~ 1 Gyr. When atmospheric CO 2 is initially high, rapid weathering draws CO 2 out of the 
atmosphere on a 1-10 Myr timescale, dramatically cooling the climate (Figures 7A &: B). This rapid 
CO 2 drawdown also ensures that the early moist greenhouse state caused by the high Pco 2 initial 
condition is short lived, and any water loss during this time would be kept to a minimum (see also 
§3.2.2). Most of the initial CO 2 drawdown takes place via supply-limited weathering (Figure 6A), 
as the extremely high initial surface temperatures push the weathering rate into the supply-limited 
regime. During this initial CO 2 drawdown phase, carbon rapidly builds up in the plate reservoir 
and is gradually subducted into the mantle (seafloor weathering also plays a role, but is much 
smaller than the continental weathering flux (Figure 6D)). After the rapid initial draw down of 
atmospheric CO 2 , the atmosphere and ocean reservoirs reach a quasi-steady-state as arc degassing 
is balanced by continental and seafloor weathering. The evolution of carbon in the atmosphere and 
hydrosphere is then controlled by the gradual loss of carbon from the plate reservoir to the mantle. 
As the mantle is regassed by subduction, the degassing flux becomes an important contributor to 
atmospheric CO 2 , and the final steady-state is reached as the degassing and arc fluxes balance the 
terrestrial weathering and seafloor weathering fluxes. Models where half of the CO 2 is initially in 
the mantle and half in the atmosphere and ocean evolve in a similar fashion. 

A different evolution occurs when all of the planet’s CO 2 initially resides in the mantle. In this 
case carbon builds up in the atmosphere and ocean through ridge degassing. A quasi-steady-state is 
quickly reached between mantle degassing and silicate weathering. Weathering also causes carbon 
to accumulate on the seafloor, and subduct into the mantle. The final steady-state is then reached 
when the mantle has degassed a sufficient supply of carbon for the deep subduction flux to balance 
the degassing flux. The timescale for reaching steady-state is thus dictated by the plate speed for 
all three initial conditions; the plate speed sets how quickly the mantle can be regassed for the 
hot initial atmosphere models, and how quickly the mantle degasses for the cold initial atmosphere 
case. 


Despite the different initial conditions leading to different temperature and plate speed evolu¬ 
tions, the same final state is reached. This occurs for two main reasons: 1) the feedbacks inherent 
in silicate weathering dominate the effect of climate on plate speed; and 2) the influence of climate 
on plate speed has a stabilizing effect, acting to drive different initial conditions to the same final 
state. At warm temperatures plate speeds are low, suppressing both mantle degassing and arc vol- 
canism and aiding in the initial drawdown of atmospheric 002- Likewise when surface temperature 
is low, a higher plate speed increases the mantle degassing rate, helping to build CO 2 back up in 
the atmosphere. The effect of variable plate speed is small, however, and the primary reason that 
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Fig. 5.— Time evolution of the mantle reservoir (a), plate reservoir (b), ocean reservoir (c), and 
atmosphere reservoir (d) for an Earth-like model starting from three different initial conditions: 
all of the carbon in the atmosphere and ocean, i.e. Rmanfl = 0 (red line), half of the carbon in 
the mantle and half in the atmosphere and ocean, i.e. Rman,o = (blue line), and all of the 

carbon in the mantle, i.e. Rmanfi = Ctot (green line). Model parameters are listed in Tables 2 &: 3. 
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Fig. 6.— Time evolution of the terrestrial weathering flux (a), degassing flux (b), subduction flux 
(c), and seafloor weathering flux (d) for the models shown in Figure 5. The arc flux, fFgub-, and 
deep subduction flux, (1 — f)Fsub, are not shown because they are both given by the subduction 
flux multiplied by a constant. 
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Fig. 7.— Time evolution of surface temperature (a), atmospheric CO 2 (b), and plate speed (c) for 
the models shown in Figure 5. 
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different initial conditions do not lead to different final states is the weathering feedback. 


The preceding discussion can be shown more generally by analyzing the steady-state solutions 
to equations (l)-(3). In particular it can be shown that, for a given total CO 2 budget (Ctot), 
land fraction, and solar luminosity, only one steady-state solution exists, and that this steady-state 
solution is stable. In steady-state, the sum of the four carbon reservoirs must equal the total CO 2 
budget, Ctot = Ratm + Roc + Rp + Rman- Using equations (1) and (2) in steady-state and Henry’s 
Law, Ctot can be written solely as a function of Roc, 


Ctot — Roc 1 1 + 


kc-^earth 


-^earthiX fland) T 


{R OC + MH2o)gmco2 
(1 - f)Vman 




fddmelt 


{v{Roc)L) 


-1 




weather 


{Roc) 


+ Fsfw{R o 


(18) 


where plate speed, v, and both weathering fluxes, Fweather and Fgf^, can be written in terms of 
the ocean carbon reservoir, rather than the atmospheric carbon reservoir, using Henry’s Law. If 
the derivative of equation (18) with respect to Roc is always positive, then there can only be one 
steady-state solution. 


Taking this derivative term by term, the first term on the right hand side of equation (18) 


gives 


d{Ratm + Ro 


Mh 2 okc-^earth 


— 1 _|_ O^c-^earth flQ") 

dRoc {Roc + MH2o)‘^mc029 

which is always positive. Taking the derivative of the second term on the right hand side of equation 
(18) yields 
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which shows that only the weathering fluxes and plate speed determine the sign of this term; as /, 
fland, fd, and dmeit are all independent of Roc, changing any of these parameters does not lead to 
multiple steady-states. Fweather increases with Roc, because a higher ocean carbon reservoir means 
higher atmospheric CO 2 concentrations, warmer temperatures, and a larger weathering flux. In 
addition, plate speed decreases with increasing Roc because higher surface temperatures lead to 
slower plates, so d{F^oather/v)/dRoc is positive due to both the temperature and CO 2 dependence 
of the weathering flux, and the temperature dependence of the plate speed (see Appendix A.l 
for the full derivation of d{F^oather/'u)/dRoc)- Outside of the supply-limited weathering regime, 
Fweather is a much stronger function of Roc than v, so the positivity of d{Fweather/v)/dRoc is 
primarily due to the temperature and Pco 2 dependence of the weathering flux. In the supply- 
limited regime, dFweather/dRoc = 0, but the inverse relationship between Roc and v ensures that 
d{Fweather/v)/dRoc remains positive. Finally, 

d{Fgfw/v) Fff^ d / F:f^MH20 \ ( koY ( aR--^ 

dRoc V*P*cl 2 ^^oc^ V V* )\P*C 02 ) \{Mh 20 + Roc)^+^ 
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which is also always positive because a > 0. Note that Fgfw has two different dependencies on 
Roc- increasing Roc increases Pco 2 therefore Fsfw directly, but also decreases the plate speed, 
which acts to slow the seafloor weathering flux. However, as shown by equation (20), it is the 
slope of Fgfyj/v that determines whether multiple steady-state solutions exist, so the plate-speed 
effect cancels out. Therefore dCtot/dRoc is positive, as Fy^oather/v and Fsfy,/v are both increasing 
functions of Roc, and only one steady-state solution exists. Even with a plate speed dependent on 
surface temperature, initial conditions do not influence the final state models evolve to (the final 
steady-state is also shown to be stable in §A.2). An important caveat to this finding is that the 
carbon cycle model used here implicitly assumes the presence of liquid water oceans and water 
vapor in the atmosphere. If initial surface temperature and pressure conditions are beyond the 
critical point for water, it is not clear that silicate weathering, and subsequent CO 2 drawdown, 
can occur; thus, whether carbon cycling can produce habitable climates starting from such high 
temperature and pressure conditions is unknown. However, as C02-rock reactions above the liquid 
water critical point are poorly constrained, such extreme conditions are beyond the scope of this 
study. 

The only way to produce multiple steady-state solutions with the present model would be for 
plate speed to have the opposite dependence on surface temperature, where higher temperatures 
lead to faster plate speeds. In this case both a high surface temperature, high plate speed with rapid 
degassing, and a low surface temperature, slow plate speed with low degassing steady-state could 
exist. However, plate speed would have to be a very strong function of surface temperature in order 
to counteract the temperature and Pco 2 dependence of the weathering fluxes and cause dCtot /dRoc 
to change sign. Furthermore, the result that only one steady-state solution exists is robust to 
uncertainties in the parameters /3, a, Fa, and a, because the sign of all these parameters is well 
known; in other words the uncertainties in these parameters are not sufficient to introduce hysteresis 
to the system because they would not change the sign of either dFy,cather/dRoc or dFgfw/dRoc- 


3.2. Transition to Supply-Limited Weathering at Small Land Areas and Large CO 2 

Inventories 

The results in §3.1 show that as long as atmospheric temperature and pressure conditions are 
below the critical point for liquid water, such that silicate weathering can occur, initial atmospheric 
CO 2 concentrations do not influence the final state reached by an Earth-sized planet with a given 
exposed land area, total planetary CO 2 budget, and solar insolation. The dependence of plate speed 
on surface temperature does not lead to multiple steady-states or hysteresis in the coupled carbon 
cycle-plate tectonic system. Furthermore, the dependence of plate speed on surface temperature 
does not diminish the ability of the carbonate-silicate cycle to buffer surface temperatures against 
changes in solar luminosity; in fact, climate stabilization is weakly enhanced because degassing 
is decreased at high surface temperatures. However, the area of exposed land and the planetary 
CO 2 budget can have a major impact on whether hospitable surface temperatures, and the climate 
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buffering effects of the long-term carbon cycle, can be maintained on a planet. 


In particular, if weathering becomes globally supply-limited a planet will not be able to regu¬ 
late atmospheric CO 2 levels through continental weathering; only the weaker seafloor weathering 
feedback would be active. In this case the planet may enter a hot, moist greenhouse climate where 
water is lost rapidly to space, either due to high atmospheric CO 2 concentrations or the inability 
of the carbon cycle to drawdown atmospheric CO 2 as solar luminosity increases. In this section, 
I show that at small land fractions and large planetary CO 2 inventories, weathering does become 
supply-limited, leading to planets with hot climates and rapid water loss; with no mechanism for 
decreasing Pco 2 these planets may lose their water and become inhospitable to life. 


The transition to a globally supply-limited weathering regime at small land fractions is shown 
in Figure 8. While weathering is kinetically-limited (i.e. for fund ~ 0.3 — 1.5 x 10“^), smaller 
land fractions lead to a modest increase in both Pco 2 surface temperature. However, once 
weathering becomes supply-limited, at land fractions lower than 1.5 x 10“^, Pc 02 surface 
temperature increase sharply as fund decreases, before flattening out as the majority of the planet’s 
total CO 2 inventory accumulates in the atmosphere and ocean. The surface temperature and 
atmospheric CO 2 trends, and the transition to supply-limited weathering, can be explained as 
follows. In the kinetically-limited regime, the change in Pco 2 ™ response to a decrease in fund is 
governed by 
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where the ridge and arc degassing fluxes can be considered constant in this regime. Therefore 
when fland decreases, Pco 2 (and in turn T) must increase in order to balance the degassing fluxes. 
However, as the terrestrial weathering flux is a strong function of Pco 2 ) only a modest increase in 
atmospheric CO 2 is needed to balance the ridge and arc degassing fluxes. The seafloor weathering 
flux also increases with increasing Pco 2 1 helping to balance the degassing fluxes with only a small 
increase in atmospheric CO 2 . However, as seafloor weathering is a weak function of Pc 02 ^ it plays 
a minor role compared to changes in the terrestrial weathering flux. 


The sharp increase in atmospheric CO 2 at fund ~ 1-5 x 10“^ represents the transition to 
globally supply-limited weathering. This transition point occurs at a critical Pco 2 > where Fweather ~ 
Pws] approximating this condition as Fweather = 0.99Tu,^, 


Pc02 

p* 

.^C02, 


exp 


Fa T idnwLwti 

Ra 


1 


= - In (0.01) 


Fmax PrXcc-'^earth fland 


F*mcc 


(23) 


which can be solved numerically to give the critical Pco 2 where weathering is supply-limited. Note 
that the critical Pco 2 is not a function of the land fraction, which cancels out in equation (23). 
When terrestrial silicate weathering is supply-limited, it is fixed at Fw^ and therefore no longer a 
function of atmospheric CO 2 . As a result, the only way to balance the degassing fluxes is through 
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Fig. 8.— Steady-state surface temperature (a) and partial pressure of atmospheric CO 2 (b) as 
a function of land fraction. The boundary between the kinetically-limited and supply-limited 
weathering regimes is shown with a dashed line, and the critical Pco 2 ^ where weathering becomes 
supply-limited, is shown in panel (b). A total CO 2 inventory of 10^^ mol is assumed, and an 
insolation of S*, the solar constant, is used. All other model parameters are listed in Tables 2 & 3. 
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a higher seafloor weathering flux and by depleting the mantle and plate reservoirs, and hence 
lowering the CO 2 degassing fluxes. The weak dependence of seafloor weathering on Pco 2 means 
that depleting the mantle and plate reservoirs, thus producing a large increase in atmospheric CO 2 , 
is the primary mechanism for bringing the system back into equilibrium. As a result, atmospheric 
CO 2 concentration increases sharply as the land fraction drops below ~ 1.5 x 10“^ (Figure 8B). A 
stronger Pc 02 dependence to seafloor weathering, possible if it is directly temperature-dependent, 
can potentially balance degassing without significant atmospheric CO 2 buildup when continental 
weathering is supply-limited. However, as shown in §3.2.3, once seafloor weathering also becomes 
supply-limited the CO 2 drawdown rate will again be unable to balance the degassing rate until the 
mantle and plate reservoirs are depleted, and hot climates will form. 

The transition between globally supply- and kinetically-limited weathering also depends on the 
planetary CO 2 inventory (Figure 9). When Ctot decreases, atmospheric CO 2 concentrations also 
decrease in response to lower mantle and arc degassing rates (owing to a smaller amount of CO 2 
throughout the mantle and surface reservoirs). As a result, the boundary between the kinetically- 
limited and supply-limited weathering regimes drops to lower fiand as Ctot decreases (Figure 9A). 
Below Ctot ~ 5 X 10^^, the supply-limited regime disappears entirely because seafloor weathering 
does not allow enough CO 2 to build up in the atmosphere for the critical Pco 2 to be reached, 
even for a planet with an infinitesimally small land fraction. Planets with large CO 2 inventories 
and small land fractions therefore lie in the supply-limited weathering regime, resulting in high 
atmospheric CO 2 concentrations and hot climates, while planets with either low CO 2 inventories or 
large land areas, like the Earth, lie in the kinetically-limited weathering regime where the carbon 
cycle can maintain a temperate climate. Overall, the kinetically-limited regime dominates, at least 
for the baseline value of Emax = 10 mm/yr chosen here, and thus the carbon cycle should be able 
to stabilize climate on most rocky planets. 


3.2.1. Influence of Erosion Rate and Sensitivity to Model Parameters 

The boundary between the supply-limited and kinetically-limited weathering regimes is strongly 
sensitive to the maximum erosion rate {Emax) used to define the supply limit to weathering, and 
less sensitive to other model parameters (Figure 10). A lower Emax means a lower supply limit 
to the weathering flux; as a result, the supply-limited weathering regime expands significantly in 
Ctot — hand space. Thus, the ability to sustain high erosion rates is important for effective cli¬ 
mate stabilization through the global carbon cycle. Given that uplift and orogeny are crucial for 
producing high erosion rates, tectonics may play a key role in keeping planets in the kinetically- 
limited weathering regime, and therefore enabling the long-term carbon cycle to maintain habitable 
climates over geologic timescales (see §4.4). 

The parameters fl and Ea, which describe the direct dependence of atmospheric CO 2 and 
surface temperature, respectively, on the rate of the weathering reaction, have a relatively weak 
influence on the boundary between the kinetically-limited and supply-limited regimes (Figures 
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Fig. 9.— Steady-state partial pressure of atmospheric COg (a) and surface temperature (b) as a 
function of land fraction and total planetary CO 2 budget. The boundary between the kinetically- 
limited and supply-limited weathering regimes is marked by a solid line, and the approximate 
position of the present day Earth is shown in panel (b). Insolation is assumed equal to the solar 
constant, and all other model parameters are listed in Tables 2 &: 3. 
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Fig. 10.— Boundary between the kinetically-limited (KL) and supply-limited (SL) weathering 
regimes for different Emax (a), j3 (b), Ea (c), / (d), dmeit (e), a (f), E*^^ (g), and S (h); for all 
plots the red curve indicates the standard parameter value used throughout the paper (e.g. in 
Figures 8 & 9). All model parameters are as specified in Tables 2 & 3, except when explicitly 
varied. The reference seafloor weathering rate, E*^^, is given in units of mol/Myr in panel (g). 
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lOB &: C). Lower values of (5, which are applicable to silicate weathering in the presence of land 
plants (e.g. Berner 1994), slightly decrease the size of the supply-limited regime, because a lower 
j3 increases the critical CO 2 level where weathering becomes supply-limited (see equation (23)), 
therefore allowing kinetically-limited weathering to prevail at smaller land fractions and larger 
planetary CO 2 inventories. Decreasing the activation energy of the weathering reaction, Ea, also 
shrinks the supply-limited regime for the same reason; the critical Pco 2 decreases at lower Ea- Both 
larger / and larger dmeit expand the supply-limited weathering regime (Figures lOD & E). When 
/, the fraction of subducted carbon that degasses at arc volcanoes, is larger, atmospheric CO 2 
levels are higher and the supply limit to weathering is reached at larger land fractions or smaller 
planetary CO 2 inventories. Likewise, a greater dmeit, the depth of melting at mid-ocean ridges, 
increases the mantle degassing rate and therefore causes higher atmospheric CO 2 concentrations. 
The influence of / and dmeit highlight the role of mantle temperature in maintaining kinetically 
controlled weathering, as both / and dmeit are in principle functions of mantle temperature (where 
higher mantle temperatures increase both / and dmeit)- 

The parameters governing seafloor weathering, a, which describes the dependence on Pco 2 , 
and the reference seafloor weathering rate, also have important influences on the weathering 

regime boundary (Figures lOF &: G). Both smaller a and smaller expand the supply-limited 
regime, because less efficient seafloor weathering means higher atmospheric CO 2 levels, and the 
critical Pco2 where weathering becomes supply-limited is reached at larger fiand or smaller Ctot- 
The case for F*j^ = 0 shows where the boundary between the weathering regimes would lie without 
any seafloor weathering. Finally, varying solar luminosity has only a very minor effect on the weath¬ 
ering regime boundary (Figure lOH); the supply-limited weathering regime is slightly expanded as 
solar luminosity increases, mainly due to a decrease in the critical Pco 2 higher luminosity. As a 
result, stellar evolution is not a major factor in determining whether a planet will be able to sustain 
kinetically-limited weathering. Other parameters, such as kc, the solubility of CO 2 in the oceans, 
and a, the exponent describing how changes in temperature influence precipitation rates, were 
found to have almost no effect on the boundary between the kinetically-limited and supply-limited 
weathering regimes. 


3.2.2. Climate Stabilization and Habitability with Supply-Limited Weathering 

In addition to producing hot climates at the present day solar luminosity, supply-limited weath¬ 
ering also eliminates the carbon cycle’s ability to regulate atmospheric CO 2 levels in response to 
changes in luminosity. This is important for planets that lie near the boundary with the kinetically- 
limited regime, where atmospheric CO 2 levels can still be moderate; these planets may have temper¬ 
ate surface temperatures at low solar luminosities, but as luminosity increases, the inability to draw 
down atmospheric CO 2 will cause temperatures to rise significantly. Likewise, planets within the 
supply-limited regime that sit near the inner edge of the habitable zone will experience much higher 
surface temperatures than those near the outer edge. Climate stabilization in the kinetically- and 
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Fig. 11.— Steady-state partial pressure of atmospheric CO 2 (a) and surface temperature (b) as a 
function of land fraction and total planetary CO 2 budget for S = 0.85"* and S = 1.25* (panels c & 
d). Model parameters are listed in Tables 2 &: 3. 
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supply-limited weathering regimes is shown in Figure 11. Within the kinetically-limited weathering 
regime, increasing solar luminosity is counteracted by a decrease in Pco 2 i thanks to the weather¬ 
ing feedback (Figures llA &: C). As a result, surface temperatures remain moderate (Figures IIB 
& D). However in the supply-limited regime, atmospheric CO 2 cannot adjust to changes in solar 
luminosity, and temperatures increase more sharply as luminosity increases. 

The hot, CO 2 rich atmospheres that result from supply-limited weathering have negative im¬ 
plications for planetary habitability. First, the typical surface temperatures estimated for the 
supply-limited regime, 400-600 K, are higher than any known life form on the modern Earth can 
tolerate (Takai et al. 2008), and thus life may not be possible at these conditions. Second, although 
these surface temperatures are within the liquid water stability held, they would place planets in 
a moist greenhouse climate state where photo dissociation of H 2 O and hydrogen escape leads to 
rapid water loss (e.g. Kasting 1988). As a result, planets in the supply-limited weathering regime 
could lose their oceans, rendering them uninhabitable. However, whether complete water loss will 
occur in practice is unclear because as sea level drops more land will be exposed, enhancing conti¬ 
nental weathering and potentially drawing enough CO 2 out of the atmosphere to stop water loss. If 
planets that originally lie in the supply-limited weathering regime can reestablish temperate, stable 
climates with higher land fractions before their oceans are completely lost to space, they will once 
again be habitable (i.e. the waterworld self-arrest mechanism proposed by Abbot et al. (2012)). 

A detailed model of water loss is beyond the scope of this study, so the likelihood of water- 
world self-arrest on supply-limited planets is assessed using a simple order of magnitude estimate. 
Specifically, the timescale for complete water loss is compared to the timescale for CO 2 drawdown 
via silicate weathering. If Earth were completely covered in oceans, the timescale for water loss 
is estimated at ~ 100 Myr (Abbot et al. 2012). To determine the timescale for CO 2 drawdown, I 
assume a planet that initially lies within the supply-limited regime, and therefore has a hot, CO 2 
rich atmosphere, undergoes water loss to a given land fraction. The weathering timescale, Tw, is 
therefore 

_ RatmjT) + RocjT) 

{l/2)F^Mand)+F,fUT)- ^ ’ 

Carbon dioxide dissolved in the oceans is included because rapid equilibration between the atmo¬ 
sphere and ocean means that CO 2 from both reservoirs must be drawn down to significantly cool 
climate. Eurthermore, continental weathering is supply-limited for all land fractions, and thus the 
CO 2 drawdown rate is primarily controlled by the land fraction, because the planets being consid¬ 
ered here have CO 2 rich atmospheres. As shown by equation (23), the transition to supply-limited 
weathering is determined solely by FCO 21 is independent of fiand- Weathering can only decrease 
from the supply-limit as atmospheric CO 2 is drawn down, so using the supply-limit to continental 
weathering in equation (24) gives the most optimistic estimate. 

Eor a planet with T = 600 K, does not reach 100 Myr until the land fraction is ~ 0.1 
(Figure 12). At this point significant water loss will have occurred, and the remaining oceans may 
be lost before weathering can cool the climate (i.e. the water loss timescale at this point will be 
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Fig. 12.— Weathering timescale, equation (24), as a function of land fraction for three different 
surface temperatures with a present day solar luminosity: T = 400 K (green line), T = 500 K (blue 
line), and T = 600 K (red line). These surface temperatures correspond to atmospheric CO 2 levels 
of Pco 2 ~ 4 bar, Pco 2 ~ 23.5 bar, and Pco 2 ~ bar, respectively. The 100 Myr timescale for 
water loss, Tiogg, is shown as a dashed line. 
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less than 100 Myrs, because most of the oceans have already been lost). The atmospheric and 
oceanic CO 2 reservoirs will also have declined by the time fiand reaches 0.1, but given the long 
weathering timescales this effect is likely negligible. However, with T = 400 K, corresponding to 
a planet that originally sat near the boundary between supply- and kinetically-limited weathering, 
the weathering timescale drops below 100 Myr at land fractions of ~ 10“^, and waterworld self¬ 
arrest appears to be quite likely. Thus the amount of CO 2 that weathering must draw down, 
which correlates to how far into the supply-limited regime a planet lies, is crucial to determining 
whether water loss can stabilize climate before the oceans are lost. A planet’s bathymetry will 
also be important. If seafloor topography is dominated by a few tall islands, with little remaining 
topography, waterworld self-arrest will be less effective because nearly all of the planet’s water 
will be lost before a significant portion of land is exposed. However, if there are large submerged 
continents surrounded by deep basins, a significant supply of water will remain even after a large 
land fraction has been exposed. Waterworld self-arrest may be viable, especially for planets that 
originally sit near the boundary with kinetically-limited weathering, but more sophisticated models 
taking into account bathymetry are needed to place better constraints on this process. 


3.2.3. Temperature-dependent Seafloor Weathering 


In this section I show that including temperature-dependent seafloor weathering and an upper 
limit on the basalt CO 2 reservoir does not significantly impact the transition to supply-limited 
weathering, or the predicted steady-state climates in fiand - Ctot space, justifying the simpler 
model used in the rest of this study where seafloor weathering depends solely on Pco 2 ■ The upper 
bound on the basalt CO 2 reservoir, Rbasmaxi i® given by calculating the number of moles of CO 2 
that are stored by complete carbonation of the seafloor to a certain depth: 

dcarb^^Earthifl fiand) P\has 


Ttba 


^bas 


(25) 


where dcarb is the depth of complete basalt carbonation, p is the density of basalt, Xbas is the 
fraction of reactable elements in basalt (CaO, MgO, and FeO), and fhbas is the average molar mass 
of CaO, MgO, and FeO (Sleep et al. 2001 ). From Sleep et al. ( 2001 ), Xbas ~ 0 - 3 , fhbas ~ 55 g 
mol~^, and a baseline value of dcarb ~ 500 m is set based on the depth that carbonation extends to 
on the modern day Earth (Alt & Teagle 1999 ). With a negligible land fraction Rbasmax ~ 4 x 10 ^^ 
mol (see Table 4). 


In order to incorporate temperature-dependent seafloor weathering with a supply-limit, as 
determined by equation (25), the carbon cycle model (equations (l)-(3)) is modified as follows: the 
plate reservoir is separated into a basalt reservoir, Rbas^ and a sedimentary reservoir, Rged- Seafloor 
weathering is assumed to supply carbon solely to the basalt reservoir, and terrestrial weathering 
supplies carbon solely to the sedimentary reservoir. The modified carbon cycle equations are 
therefore 

^Rsed Rweather Rse dvL 
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for Rhas = Rbasmax- ^ seafloor weathering activation energy of E^f = 50 kJ mol ^ is assumed, in 
the middle of the range of activation energies given by Brady &: Gfslason (1997). 


Solving for the steady-state solutions to equations (26)-(29) as a function of fiand and Ctot 
shows that adding a direct temperature feedback and a limit to the size of the basalt reservoir does 
not significantly change the predicted climate states (Figure 13A,B). Supply-limited weathering 
still leads to extremely hot climates at low land fractions and high CO 2 budgets, because the 
basalt reservoir is not large enough to prevent the accumulation of a CO 2 rich atmosphere at these 
conditions. The supply-limited regime is slightly smaller when compared to the models without 
temperature-dependent seafloor weathering (Figure 9), with higher total CO 2 budgets required 
for weathering to become supply-limited. The difference is minor because even without a direct 
temperature feedback, seafloor weathering still acts as a major carbon sink at low land fractions and 
large CO 2 inventories. The only way to shrink the supply-limited regime further through seafloor 
weathering is with a larger carbonation depth, as shown by Figure 13C (increasing either Egf or a 
does not effect the weathering regime boundary because the boundary is dictated by the maximum 
size of the basalt CO 2 reservoir). However, dcarb > 6 km would be required to completely eliminate 
the supply-limited weathering regime for the range of Ctot and fiand explored here. Complete 
carbonation of basalt to over 6 km depth, or even beyond the top few hundred meters, is not 
likely because hydrothermal carbonate formation is generally confined to the near surface layer of 
pillow basalts where permeability is high and significant hydrothermal circulation can occur (Alt 
&: Teagle 1999). High surface temperatures and atmospheric CO 2 concentrations will make the 
accessible basalt more reactive, but are unlikely to cause deeper hydrothermal circulation because 
the seafloor permeability structure would be unaffected. Moreover, the supply-limited weathering 
regime is still strongly influenced by the maximum erosion rate, with lower erosion rates leading to 
a larger supply-limited regime, even for the modified seafloor weathering flux used in this section 
(Figure 13D). 
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Fig. 13.— Steady-state partial pressure of atmospheric CO 2 (a), and surface temperature (b) for the 
modified carbon cycle model including temperature-dependent seafloor weathering (equations (26)- 
(31); see Table 4 for parameter values). The dashed line shows when seafloor weathering becomes 
supply-limited (i.e. when R^as = Rbasmax)^ the solid line shows when continental weathering 
becomes snpply limited. Weathering is considered to be globally supply-limited when both seafloor 
and terrestrial weathering are supply-limited. The boundary between globally supply-limited (SL) 
and kinetically-limited (KL) weathering is plotted for varying dcarb (c) and Emax (d). 
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4. Discussion 

4.1. Implications for Planetary Habitability 

Overall, the results of this study have generally positive implications for planetary habitability. 
The results show that the inverse dependence of plate speed on surface temperature proposed by 
Foley & Bercovici (2014) and Foley et al. (2014) does not inhibit the carbon cycle’s ability to 
regulate climate (§3.1.1). Furthermore, the results also show that as long as silicate weathering 
can occur, initial atmospheric CO 2 concentrations, and therefore temperatures, do not influence the 
final steady-state reached (§3.1.2); planets with a given exposed land area, total CO 2 inventory, and 
receiving a given solar irradiance will evolve to the same surface temperature and Pco 2 regardless 
of the initial distribution of carbon between mantle and surface reservoirs. Thus, an initially hot 
climate is not an obstacle for developing temperate surface temperatures and active plate tectonics 
on a planet. Other factors, however, can be an obstacle for the development of a habitable climate. 
Planets with small land fractions and large CO 2 inventories are found to enter a supply-limited 
weathering regime, where weathering depends only on the erosion rate and exposed land area, and 
not on surface temperature and atmospheric CO 2 (§3.2). As a result, planets in the supply-limited 
weathering regime tend to develop hot, CO 2 rich climates because the carbonate-silicate cycle 
can no longer regulate atmospheric CO 2 levels. Such hot climates are less favorable for life and 
could lead to signihcant water loss due to photodissociation of H 2 O in the upper atmosphere and 
hydrogen escape to space. Furthermore, carbon cycling in the supply-limited weathering regime 
can no longer regulate atmospheric CO 2 levels in response to changes in solar luminosity (§3.2.2). 

Nevertheless, the supply-limited regime only takes up a relatively small portion of the fiand- 
Ctot parameter space. Thus, only planets that have significantly larger amounts of both water 
and carbon dioxide than the Earth would be at risk of losing the carbon cycle’s stabilizing climate 
feedbacks due to supply-limited weathering. As long as erosion rates are high, planets with a 
CO 2 budget similar to Earth’s can develop habitable climates, with kinetically-limited weathering, 
at land fractions as small as 10“^ — 10“^. As a result, only a handful of Hawaiian type island 
chains would be needed to maintain an active carbon cycle with stabilizing climate feedbacks (the 
Hawaiian islands alone have a land fraction of ~ 5 x 10“^). A significant percentage of Earth-like 
planets within the habitable zone should therefore be able to develop temperate climates and plate 
tectonics, and maintain these habitable surface temperatures over geologic timescales. Moreover, 
some planets that fall into the supply-limited weathering regime may even be able to reestablish 
stable, temperate climates as water loss exposes more land (§3.2.2). In particular, planets that 
originally sit near the boundary with kinetically-limited weathering only need to expose modest 
amounts of land before the timescale to cool the climate through weathering becomes shorter than 
the timescale for complete water loss. This waterworld self-arrest mechanism may not work for all 
planets, however. Planets that initially lie deep within the supply-limited regime have much higher 
atmospheric CO 2 concentrations, resulting in much longer climate cooling timescales; thus complete 
water loss may occur before silicate weathering can significantly cool the climate on these planets. 


Even so, the possibility of waterworld self-arrest further increases the odds that a significant fraction 
of planets forming within the habitable zone will indeed be habitable over a large portion of their 
geologic lifetimes. 


4 . 1 . 1 . Influence of Planet Size 

Although this study focuses solely on Earth-sized planets, some preliminary conclusions about 
rocky planets of different sizes can be drawn. Planet size influences whether plate tectonics can 
occur on a planet (e.g. O’Neill et al. 2007; Foley et al. 2012), the structure of the atmosphere and 
hence the strength of the radiative forcing provided by greenhouses gasses such as CO 2 and H 2 O 
(e.g. Kasting et al. 1993; Kopparapu et al. 2014), and the rate at which water is lost to space 
(e.g. Hunten 1973; Melosh &: Vickery 1989). The influence of size on the propensity for plate 
tectonics on a planet is controversial (Valencia et al. 2007; O’Neill & Lenardic 2007; Kite et al. 
2009; Valencia & O’Connell 2009; Korenaga 2010; van Heck & Tackley 2011; Foley et al. 2012; 
Stein et al. 2013; Noack &: Breuer 2014), because the mechanism that leads to plate tectonics on 
Earth is unknown (e.g. Tackley 2000; Bercovici 2003), and because key material properties, such 
as viscosity and thermal conductivity, are not well constrained at the temperature and pressure 
conditions expected for super-Earth planets (Karato 2011; Stamenkovic et al. 2011). 

With the grain-damage mechanism for generating plate tectonics from mantle convection, Foley 
et al. (2012) found that larger planets are more likely to have plate tectonics because larger planets 
have more vigorously convecting mantles and therefore larger driving forces for grainsize reduction. 
As a result, small planets, like Mars, may not able to develop plate tectonics even with cool surface 
temperatures. Furthermore, small planets lose their water to space more easily because of their 
smaller gravity (e.g. Kasting Sz Catling 2003). Thus small rocky planets are unlikely to maintain 
liquid water oceans because of both a lack of plate tectonics and an inability to prevent rapid water 
loss to space. On the other hand larger planets, in addition to being more likely to have plate 
tectonics, also have a larger habitable zone, because higher planet mass reduces the greenhouse 
effect of H 2 O, thereby moving the inner edge of the habitable zone to lower orbital distances 
(Kopparapu et al. 2014). However, the volatile content of super-Earths will also be important; 
if these planets are volatile rich they could end up in a supply-limited weathering regime, due 
to a small land area and a high planetary CO 2 budget. Furthermore, the conclusion from Foley 
et al. (2012) that plate tectonics is more likely on larger planets should be considered preliminary 
until better constraints are placed on the material properties of super-Earth mantles, and these 
constraints are included in numerical models of mantle convection with grain-damage. 
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4.2. Model Uncertainties 
4-2.1. Parameter Uncertainties 

The carbon cycle model used in this study relies on a large number of parameters (Table 2), 
many of which are not well constrained. However, as seen in §3, the overall results of this study 
hold for a wide range of model parameters, and are therefore robust. The result from §3.1, that the 
dependence of plate speed on surface temperature does not significantly alter the negative climate 
feedbacks of the carbon cycle, nor does it add hysteresis to the system, holds for a wide parameter 
range. The reason the dependence of plate speed on surface temperature has such a minor impact 
on the long-term carbon cycle’s climate buffering capacity is that the influence of temperature on 
the weathering flux is significantly stronger than its influence on plate speed. Either the activation 
energy for the weathering reaction {Ea) would have to be significantly lower than what is inferred 
from laboratory results and field studies, or the dependence of plate speed on surface temperature 
would have to be far stronger than what is inferred from numerical convection models, for variations 
in plate speed to significantly alter the carbon cycle. Even then, the influence of surface temperature 
on plate tectonics has a stabilizing effect, because higher surface temperatures lead to slower plates, 
and therefore lower degassing rates. Thus, some form of climate stabilization would be preserved 
as long as plate speed depends inversely on surface temperature, as expected by mantle convection 
studies with grain-damage (Eoley &: Bercovici 2014; Foley et al. 2014). 

The result of §3.1.2, that there is only one steady-state for a planet with a given land frac¬ 
tion and total CO 2 budget to evolve to, is similarly robust. This result is independent of model 
parameters, as long as the two most fundamental aspects of the model, that the weathering flux 
never decreases with increasing atmospheric CO 2 and that the plate speed decreases with increasing 
surface temperature, hold true. Note that the results of §3.1.2 still hold even when temperature- 
dependent seafloor weathering and a limit to basalt carbonation are included because the quotient of 
the seafloor weathering flux and the plate velocity would still never be a decreasing function of Roc- 
The only plausible way to introduce hysteresis into the coupled carbon cycle-plate tectonic system 
would be through additional feedbacks, not included in the present model; furthermore any new 
feedbacks would have to be very strong in order to counteract the weathering feedbacks. Finally, the 
position of the boundary between the kinetically-limited weathering and supply-limited weathering 
regimes changes with key model parameters (§3.2.1) or with the inclusion of temperature-dependent 
seafloor weathering (§3.2.3), but the overall trend of supply-limited weathering occurring at small 
land fractions and large CO 2 inventories holds for a wide range of parameters. 


4 . 2 . 2 . Uncertainties in Model Formulation 

One important caveat to the result that there is only one steady-state solution that planets 
evolve to from §3.1.2 is that a constant albedo was used in this study. Introducing a variable 
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albedo, capable of capturing the ice-albedo feedback, could lead to a situation where a planet 
oscillates back and forth between ice covered and ice free states. Such oscillations can occur during 
transient periods of anomalously high weathering rates (Mills et al. 2011), when degassing rates 
are low (Kadoya &: Tajika 2014), or when luminosity is low (Menou 2015). However, even in these 
cases hysteresis caused by different initial atmospheric CO 2 concentrations is not expected, because 
the ice free-ice covered oscillations are caused when the temperature predicted by the steady-state 
carbon cycle is lower than a threshold value, where runaway glaciation occurs. Thus, the presence 
of an oscillatory state does not affect the way a planet evolves into that state from different initial 
atmospheric CO 2 conditions, and the arguments from §3.1.2 that hysteresis does not occur still 
hold. 

Another important source of uncertainty is the climate parameterization, which is a simple 
approximation of the results of more sophisticated radiative-convective atmosphere models. A 
more advanced climate model would likely change the details of the results of this study, both 
in terms of the exact temperatures and atmospheric CO 2 concentrations calculated in §3, and 
the exact position of the weathering regime boundary in fiand — Ctot space. In particular, the 
climate parameterization used in this study likely over-predicts surface temperature at very high 
atmospheric CO 2 levels (e.g. above ~ 10 bar). Thus the exact temperatures for planets in the 
supply-limited weathering regime, or the initial temperatures predicted for planets with CO 2 rich 
primordial atmospheres, may be different from what is shown here. However, these differences 
do not impact the key results in any major way: the ability of the long-term carbon cycle to 
stabilize climate is assessed using models at much lower temperatures, where equation (15) is more 
accurate, and the transition between supply-limited and kinetically-limited weathering also occurs 
at relatively low temperatures. Finally, the lack of hysteresis caused by different initial atmospheric 
CO 2 conditions does not depend on the details of the climate parameterization. 


4.3. Comparison to Previous Studies on Weathering with Small Land Areas 

The topic of silicate weathering and climate stabilization on ocean dominated planets is rela¬ 
tively new. Abbot et al. (2012) studied this topic using a similar carbon cycle model, albeit without 
a treatment for supply-limited weathering, and found that a habitable, stable climate can still be 
maintained on an Earth-sized planet with an Earth-like CO 2 budget at land fractions as low as 
10“^. My results show that this main conclusion of Abbot et al. (2012) holds even when supply- 
limited weathering is included, as long erosion rates are higher than ~ 0.1 mm yr“^ (see Figure 
10 A). Furthermore, my results extend the study of climate stabilization via the long-term carbon 
cycle to lower land fractions and varying planetary CO 2 inventories, finding that land fractions 
must be very low, or CO 2 budgets much larger than Earth’s, for the carbon cycle to be unable to 
stabilize planetary climate. My results also provides some additional constraints on the possibil¬ 
ity of waterworld self-arrest, as proposed by Abbot et al. (2012), by showing how the weathering 
timescale varies with both land fraction and atmospheric CO 2 content. The case with T = 400 
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K (see Figure 12) is most comparable to the situation considered by Abbot et al. (2012) as the 
atmospheric CO 2 levels are similar, and shows a consistent weathering timescale of 1-10 Myr for a 
land fraction greater than 10“^. Thus supply-limited weathering does not significantly impact the 
weathering timescale, again as long as erosion rates are high. Atmospheric CO 2 content, however, 
is important as the longer weathering timescales that result from higher surface temperatures, and 
thus higher Pco 2 ^ attest to. 


4.4. Importance of Plate Tectonics for Silicate Weathering and Planetary 

Habitability 

Plate tectonics has long been considered to be an important, and perhaps even necessary, fea¬ 
ture of a habitable planet, because plate tectonics facilitates the long-term carbon cycle (e.g. Walker 
et al. 1981; Gonzalez et al. 2001; Kasting & Catling 2003). Plate tectonics rejuvenates the supply of 
weatherable rock at the surface through uplift and orogeny, and seafloor production at mid-ocean 
ridges. Plate tectonics also allows for the recycling of carbon into the mantle, which sustains CO 2 
degassing to the atmosphere at both arcs and ridges. By including the effects of supply-limited 
weathering, this study directly illustrates the importance of plate tectonics in planetary habitability, 
through the influence of plate tectonics on erosion rates. The maximum erosion rate, which controls 
the supply limit to weathering, is related to tectonics because over long time scales erosion can only 
proceed as quickly as tectonics brings rock to the surface. Uplift and volcanic resurfacing therefore 
provide the ultimate, long-term upper bound on erosion rates. As orogeny, uplift, and volcanism 
are essential aspects of plate tectonics, planets with plate tectonics should have higher maximum 
erosion rates, and therefore smaller supply-limited weathering regimes. In other words, planets 
with plate tectonics are more likely to have a carbon cycle that regulates climate, and are therefore 
more likely to be habitable, because higher erosion rates mean kinetically-limited weathering can 
be sustained at smaller land fractions and larger planetary CO 2 inventories. Furthermore, as plate 
tectonics is primarily responsible for the formation of continental crust on Earth via arc volcanism, 
plate tectonic planets may also have higher land fractions compared to stagnant lid planets. 


5. Conclusions 

A simple coupled plate tectonic-carbon cycle model shows that the ability of the long-term 
carbon cycle to buffer climate in response to changes in solar luminosity is preserved when the 
dependence of plate tectonics on surface temperature, based on recent geodynamical studies, is 
included. The climate feedbacks involved in silicate weathering far outweigh the influence climate 
has on plate tectonics. Moreover, the influence of surface temperature on plate speed actually 
enhances climate buffering, because warmer temperatures lead to slower plate speeds and lower 
CO 2 degassing rates, in addition to higher weathering rates. Furthermore, initial conditions are 
found to not influence the final state reached when the carbon cycle comes to steady-state. As long 
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as liquid water is present and silicate weathering can occur, a planet with a given exposed land area, 
total CO 2 budget, and receiving a given solar irradiance will evolve to the same surface temperature, 
regardless of initial atmospheric CO 2 concentrations. An initially hot, CO 2 rich atmosphere is not 
an obstacle for developing a temperate climate and active plate tectonics on a planet. 

However, the area of exposed land and total CO 2 budget of a planet both influence whether 
the steady-state climate that planet evolves to is clement. At low land areas and high CO 2 bud¬ 
gets continental weathering becomes supply-limited, meaning the supply of bedrock to the surface 
via erosion controls the weathering rate, rather than the kinetics of the reaction between CO 2 
and silicate minerals. When weathering is supply-limited, it no longer depends on atmospheric 
CO 2 concentration, and the climate feedback is lost. Thus, planets that lie in the supply-limited 
weathering regime have hot, CO 2 rich climates and lack the weathering feedbacks necessary to 
draw down atmospheric CO 2 concentrations in response to increases in solar luminosity. The 
hot climates that result from supply-limited weathering will cause rapid water loss. Whether the 
oceans are completely lost, however, likely depends on the amount of CO 2 in the atmosphere of a 
planet in the supply-limited regime. Planets lying near the edge of the supply-limited regime have 
lower atmospheric CO 2 concentrations and can probably reestablish temperate, stable climates as 
water loss exposes more land and enhances terrestrial weathering. However, planets lying deep 
within the supply-limited regime have much larger atmospheric CO 2 concentrations and therefore 
may undergo total water loss before weathering can re-stabilize the climate. Thus ocean domi¬ 
nated planets with CO 2 inventories significantly greater than Earth’s, such that they would lie well 
within the supply-limited regime, are less likely to be habitable. High erosion rates are important 
for preventing supply-limited weathering on a planet because they enhance the supply of bedrock 
to the surface. The role of erosion rate further highlights the importance of plate tectonics for 
planetary habitability: plate tectonics drives uplift and orogeny, and therefore plays a direct role 
in maintaining high erosion rates on a planet. 


I thank Josh West for discussions on supply-limited weathering, John D. Platt for help with 
mathematical analysis of the steady-state solutions, and John Chambers and Peter Driscoll for 
comments on an earlier version of the manuscript. I also thank Dorian Abbot for a thoughtful 
and constructive review that helped to improve the manuscript. This work was supported by the 
NASA Astrobiology Institute under cooperative agreement NNA09DA81A. 
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Table 4; Table of variables and parameters for temperature-dependent seafloor weathering. 


Symbol 

Definition 

Baseline Value 

Equation 

dcarb 

Depth of complete basalt carbonation 

500 m (SOI) 

(25) 

p 

Basalt density 

2800 kg m-3 (TS) 

(25) 

Xbas 

Fraction of CaO, MgO, and FeO in basalt 

0.3 (SOI) 

(25) 

'^bas 

Molar mass of CaO, MgO, and FeO 

55 g mol~^ (derived) 

(25) 

d^basmax 

Maximum size of basalt reservoir 

4 X 10^^ mol (SOI) 

(25) 

Rsed 

Seafloor sedimentary carbon reservoir 

- 

(26) 

Rbas 

Seafloor basalt carbon reservoir 

- 

(27) 

FsfWT 

Temperature-dependent seafloor weathering flux 

- 

(30) 

Esf 

Activation energy for seafloor weathering 

50 kJ mol“^ (B97) 

(30) 


Key for citations: B97=(Brady &: Gi'slason 1997), S01=(Sleep et al. 2001). 
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A. Appendix 

A.l. Weathering Flux Derivative 

The derivative of F^eather/v with respect to Roc, as introduced in equation (20), is derived in 
this section. Using equations (11) and (14), 
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and Vs = 0.011 cm yr ^ K ^ (see equation (17)). Evaluating the second term on the right hand 
side of equation (Al), 
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Combining equations (Al), (A2), and (A3) shows that the derivative of Fweather/v with respect to 
Roc is positive, as described in the main text (§3.1.2). 


A.2. Linear Stability Analysis 


The stability of the steady-state solution for a given planetary CO 2 budget, land fraction, 
and solar luminosity can be assessed using a linear stability analysis. Although the global carbon 
cycle model involves three differential equations, there are only two independent equations because 
one reservoir can always be related to the other reservoirs via the planetary CO 2 budget: Ctot = 
Ratm+Roc + Rp + Rman- I chose to treat Rp and Rman as the independent dimensions, and calculate 
the atmosphere and ocean reservoirs from Ratm + Roc = Ctot — Rp — Rman- The stability of the 
resulting two dimensional system of ordinary differential equations is tested using the Jacobian 
matrix, J, 
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where the derivatives of / and g are evaluated at Rp and Rman, the steady-state values of the plate 
and mantle reservoirs, respectively (hats denote a quantity at steady-state throughout this section). 
If the trace of J is negative and the determinant is positive, then the steady-state solution is stable 
(Strogatz 1994). All terms in / and g that are functions of Pco 2 ^ effectively functions 
of Rp and Rman through the planetary CO 2 budget; i.e. perturbations in Rp and Rman also give 
rise to perturbations in Ratm and Roc in order to keep the planet’s CO 2 budget balanced. The 
atmospheric carbon reservoir is written in terms of Rp and Rman as Ratm + R oc — Ctot Rp Rman-) 
where 

Ratm'^C029 _ kcRoc (A7) 

-^earth Roc + MH20 

from Henry’s law. The amount of CO 2 dissolved in the ocean is typically small compared to the 
number of moles of H 2 O, meaning Roc « Mh 20 'i using this approximation 
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because Fs^pl — /) = Fdegas at steady-state, and 
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Therefore df /dRman is guaranteed to be negative, and df /dRp is positive. The derivatives of g 
are more complicated due to the seafloor and terrestrial weathering fluxes, but can be written as 
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Atmospheric CO 2 decreases with increasing Rp or Rman, because increasing the amount of carbon 
in the plate or mantle reservoirs leads directly to a decrease in the atmospheric and oceanic carbon 
reservoirs in order to keep the CO 2 budget balanced. As a result, temperature also decreases with 
increasing Rp or Rman, and plate speed increases due to the drop in surface temperature. As chang¬ 
ing either the plate or mantle carbon reservoirs causes a corresponding change in the atmospheric 
reservoir, Ratm is a linear function of both Rp and Rman with the same slope (see equation (A8)), 
and dPco2/dRp = dPco2/dRman, dT/dRp = dT/dRman, and dv/dRp = dv/dRman- 


Using equations (A13)-(A15), the derivatives of the terrestrial and seafloor weathering fluxes 
can be determined. As derivatives of Pco 2 , T, and v with respect to Rp are equal to those with 
respect to Rman, the derivatives of Fweather and Fgfw with respect to Rp are also equal to those 
with respect to Rman- Evaluating those derivatives gives 
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The fact that plate speed increases with increasing Rp or Rman ensures that the first term on 
the right hand side of both equations (All) and (A12) is negative; a positive perturbation in the 
plate reservoir will increase the subduction flux both because the plate reservoir is larger, and the 
plate speed is higher, thereby acting to deplete the plate reservoir and damp out the perturbation. 
Likewise, a positive perturbation in the mantle reservoir also increases the plate speed and hence 
acts to decrease the plate reservoir. The terrestrial weathering flux decreases with increasing Rp or 
Rman (see equation (A17)) because both Pco 2 T are lowered by an increase in either the plate 
or mantle carbon reservoirs. As a result, the second term on the right hand side of equations (All) 
and (A12) is also negative. There are two competing effects that determine the sign of dFgjru/dRp 
(see equation (A16)); increasing Rp drops Pc 02 , acting to decrease Fgfrv, but also increases v, acting 
to increase Fgfru- The influence of decreasing atmospheric CO 2 tends to dominate, so dFgj^/dRp 
is found to be negative for all parameter ranges explored in this study. 
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Fig. Al.— Logarithm of the negative of the trace of the Jacobian, J, (a) and logarithm of the 
determinant of J (b) for steady-state solutions to equations (l)-(3) as a function of Ctot and fiand- 
A present day solar luminosity is assumed, and all other parameters are listed in Tables 2 & 3. 
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Summing all of the terms, dg/dRp and dg/dRman are negative. Thus the trace of J, d f / dRman+ 
dg/dRp, is negative, and the determinant, {d f/dRman){dg/ dRp) — {dg /dRman){df /dRp) , is pos¬ 
itive, and steady-state solutions are stable. Figure A1 shows the logarithm of the negative of the 
trace of J and the logarithm of the determinant of J for steady-state solutions to equations (l)-(3) 
over a wide range of Ctot and fiand- Though the trace and determinant approach 0 at large Ctot 
and small fiand^ particularly in the supply-limited weathering regime (see §3.2), the trace remains 
negative and determinant remains positive, demonstrating that these steady-states are stable. 


-48- 


REFERENCES 

Abbot, D. S., Cowan, N. B., &: Ciesla, F. J. (2012). Indication of Insensitivity of Planetary 
Weathering Behavior and Habitable Zone to Surface Land Fraction. Astrophys. J., 756, 
178. 

Abbot, D. S., Voigt, A., & Roll, D. (2011). The Jormungand global climate state and implications 
for Neoproterozoic glaciations. J. Geophys. Res. Atmospheres, 116, 18103. 

Ague, J. J., & Nicolescu, S. (2014). Carbon dioxide released from subduction zones by fluid- 
mediated reactions. Nature Geoscience, 7, 355-360. 

Alt, J. C., & Teagle, D. A. H. (1999). The uptake of carbon during alteration of ocean crust. 
Geoehim. Cosmoehim. Aeta, 63, 1527-1535. 

Bercovici, D. (2003). The generation of plate tectonics from mantle convection. Earth Planet. Sci. 
Lett, 205, 107-121. 

Bercovici, D., & Ricard, Y. (2014). Plate tectonics, damage, and inheritance. Nature, 508, 513-516. 

Berner, R. (2004). The Phanerozoic carbon cycle. Oxford University Press. 

Berner, R. A. (1994). GEOCARB II; a revised model of atmospheric CO 2 over Phanerozoic time. 
Am. J. Sei., 294 , 56-91. 

Berner, R. A., Lasaga, A. C., & Carrels, R. M. (1983). The carbonate-silicate geochemical cycle 
and its effect on atmospheric carbon dioxide over the past 100 million years. Am. J. Sei, 
283, 641-683. 

Brady, P. V. (1991). The effect of silicate weathering on global temperature and atmospheric CO 2 . 
J. Geophys. Res., 96, 18101. 

Brady, P. V., &: Gi'slason, S. R. (1997). Seafloor weathering controls on atmospheric CO 2 and 
global climate. Geoehim. Gosmochim. Acta, 61, 965-973. 

Catling, D. C., & Claire, M. W. (2005). How Earth’s atmosphere evolved to an oxic state: A status 
report. Earth Planet. Sci. Lett., 237, 1-20. 

Charnay, B., Forget, F., Wordsworth, R., Leconte, J., Millour, E., Codron, F., &: Spiga, A. (2013). 
Exploring the faint young Sun problem and the possible climates of the Archean Earth with 
a 3-D GCM. J. Geophys. Res. Atmospheres, 118, 10414. 

Christensen, U. (1989). Mixing by time-dependent convection. Earth Planet. Sci. Lett., 95, 382- 
394. 

Coogan, L. A., &: Dosso, S. E. (2015). Alteration of ocean crust provides a strong temperature 
dependent feedback on the geological carbon cycle and is a primary driver of the sr-isotopic 
composition of seawater. Earth Planet. Sci. Lett., 4^5, 38-46. 



- 49 - 


Coogan, L. A., &: Gillis, K. M. (2013). Evidence that low-temperature oceanic hydrothermal 
systems play an important role in the silicate-carbonate weathering cycle and long-term 
climate regulation. Geochem. Geophys. Geosyst., I 4 , 1771-1786. 

Dasgupta, R., &: Hirschmann, M. M. (2010). The deep carbon cycle and melting in Earth’s interior. 
Earth Planet. Sci. Lett., 298, 1-13. 

de Pater, I., & Lissauer, J. (2010). Planetary Sciences. (2nd ed.). New York: Cambridge Univ. 
Press. 

Delano, J. W. (2001). Redox History of the Earth’s Interior since ~3900 Ma; Implications for 
Prebiotic Molecules. Origins Life Evol. B., 31, 311-341. 

Driscoll, P., & Bercovici, D. (2013). Divergent evolution of earth and venus; Influence of degassing, 
tectonics, and magnetic fields. Icarus, 226, 1447-1464. 

Edmond, J. M., Palmer, M. R., Measures, C. I., Grant, B., & Stallard, R. F. (1995). The fluvial 
geochemistry and denudation rate of the Guayana Shield in Venezuela, Colombia, and Brazil. 
Geochim. Cosmochim. Acta, 59, 3301-3325. 

Feulner, G. (2012). The faint young Sun problem. Rev. Geophys., 50, 2006. 

Foley, B. J., & Bercovici, D. (2014). Scaling laws for convection with temperature-dependent 
viscosity and grain-damage. Geophys. J. Int., 199, 580-603. 

Foley, B. J., Bercovici, D., &: Elkins-Tanton, L. T. (2014). Initiation of plate tectonics from post¬ 
magma ocean thermochemical convection. ,J. Geophys. Res. Solid Earth, 119, 8538-8561. 

Foley, B. J., Bercovici, D., & Landuyt, W. (2012). The conditions for plate tectonics on super¬ 
earths: Inferences from convection models with damage. Earth Planet. Sci. Lett., 331332, 
281 - 290. 

Fowler, C. M. R. (2005). The Solid Earth: An Introduction to Global Geophysics. (2nd ed.). New 
York: Cambridge Univ. Press. 

Franck, S., Block, A., von Bloh, W., Bounama, C., Schellnhuber, H.-J., & Svirezhev, Y. (2000). 
Habitable zone for Earth-like planets in the solar system. Planet. Space Sci., 48 , 1099-1105. 

Franck, S., Kossacki, K., & Bounama, C. (1999). Modelling the global carbon cycle for the past 
and future evolution of the earth system. Ghem. GeoL, 159, 305-317. 

Frost, D. J., Mann, U., Asahara, Y., &: Rubie, D. C. (2008). The redox state of the mantle during 
and just after core formation. Phil. Trans. R. Soc. A, 366, 4315-4337. 

Frost, D. J., &: McCammon, C. A. (2008). The Redox State of Earth’s Mantle. Annu. Rev. Earth 
Planet. Sci., 36, 389-420. 



- 50 - 


Gabet, E. J., & Mudd, S. M. (2009). A theoretical model coupling chemical weathering rates with 
denudation rates. Geology, 37, 151-154. 

Gaillardet, J., Dupre, B., Louvat, P., &: Allegre, C. (1999). Global silicate weathering and co 2 
consumption rates deduced from the chemistry of large rivers. Chem. GeoL, 159, 3-30. 

Gillis, K. M., &: Googan, L. A. (2011). Secular variation in carbon uptake into the ocean crust. 
Earth Planet. Sei. Lett., 302, 385-392. 

Gonzalez, G., Brownlee, D., & Ward, P. (2001). The Galactic Habitable Zone: Galactic Chemical 
Evolution. learus, 152, 185-200. 

Gough, D. O. (1981). Solar interior structure and luminosity variations. Solar Phys., 74, 21-34. 

Halliday, A. N., & Wood, B. J. (2007). The composition and major reservoirs of the Earth around 
the time of the moon-forming impact. In D. Stevenson, & G. Schubert (chief editor) (Eds.), 
Treatise On Geophysics (pp. 13-50). New York: Elsevier volume 9. 

Haqq-Misra, J. D., Domagal-Goldman, S. D., Kasting, P. J., & Kasting, J. F. (2008). A Revised, 
Hazy Methane Greenhouse for the Archean Earth. Astrobiology, 8, 1127-1137. 

Hart, M. H. (1978). The evolution of the atmosphere of the earth. Icarus, 33, 23-39. 

Hart, M. H. (1979). Habitable zones about main sequence stars. Icarus, 37, 351-357. 

Hilley, G. E., Chamberlain, C. P., Moon, S., Porder, S., & Willett, S. D. (2010). Competition 
between erosion and reaction kinetics in controlling silicate-weathering rates. Earth Planet. 
Sei. Lett, 293, 191-199. 

Hoffman, N. R. A., & McKenzie, D. P. (1985). The destruction of geochemical heterogeneities by 
differential fluid motions during mantle convection. Geophys. J. Roy. Astr. S., 82, 163-206. 

Hunten, D. M. (1973). The Escape of Light Gases from Planetary Atmospheres. J. Atmos. Sei., 
30, 1481-1494. 

Kadoya, S., & Tajika, E. (2014). Conditions for Oceans on Earth-like Planets Orbiting within the 
Habitable Zone: Importance of Volcanic CO 2 Degassing. Astrophys. J., 790, 107. 

Karato, S. (2011). Rheological structure of the mantle of a super-earth: Some insights from mineral 
physics. Icarus, 212, 14 - 23. 

Kasting, J. F. (1988). Runaway and moist greenhouse atmospheres and the evolution of earth and 
Venus. Icarus, 74, 472-494. 

Kasting, J. F. (1989). Long-term stability of the earth’s climate. Glob. Planet. Ghange, 1, 83-95. 

Kasting, J. F., & Ackerman, T. P. (1986). Climatic consequences of very high carbon dioxide levels 
in the earth’s early atmosphere. Science, 234, 1383-1385. 



- 51 - 


Kasting, J. F., & Catling, D. (2003). Evolution of a Habitable Planet. Annu. Rev. Astron. Astro- 
phys., 41 , 429-463. 

Kasting, J. F., Pollack, J. B., & Ackerman, T. P. (1984). Response of earth’s atmosphere to 
increases in solar flux and implications for loss of water from Venus. Icarus, 57, 335-355. 

Kasting, J. F., Whitmire, D. P., & Reynolds, R. T. (1993). Habitable Zones around Main Sequence 
Stars. Icarus, 101, 108-128. 

Kelemen, P. B., &: Manning, C. E. (2015). Reevaluating carbon fluxes in subduction zones, what 
goes down, mostly comes up. Proc. Natl Acad. Sci., 112, E3997-E4006. 

Kellogg, L. H., &: Turcotte, D. L. (1990). Mixing and the distribution of heterogeneities in a 
chaotically convecting mantle. J. Geophys. Res., 95, 421-432. 

Kite, E. S., Manga, M., & Gaidos, E. (2009). Geodynamics and Rate of Volcanism on Massive 
Earth-like Planets. Astrophys. J., 700, 1732-1749. 

Kleine, T., Miinker, G., Mezger, K., & Palme, H. (2002). Rapid accretion and early core formation 
on asteroids and the terrestrial planets from Hf-W chronometry. Nature, 4^^: 952-955. 

Kopparapu, R. K., Ramirez, R. M., SchottelKotte, J., Kasting, J. P., Domagal-Goldman, S., & 
Eymet, V. (2014). Habitable Zones around Main-sequence Stars: Dependence on Planetary 
Mass. Astrophys. J., 787, L29. 

Korenaga, J. (2006). Archean geodynamics and the thermal evolution of earth. In K. Benn, J.-G. 
Mareschal, & K. Condie (Eds.), Archean Geodynamics and Environments (pp. 7-32). AGU 
Geophysical Monograph Series volume 164. 

Korenaga, J. (2010). On the likelihood of plate tectonics on super-earths: Does size matter? 
Astrophys. J., 725, L43. 

Kump, L. R., Brantley, S. L., & Arthur, M. A. (2000). Ghemical Weathering, Atmospheric GO 2 , 
and Glimate. Annual Review of Earth and Planetary Sciences, 28, 611-667. 

Landuyt, W., &: Bercovici, D. (2009). Variations in planetary convection via the effect of climate 
on damage. Earth Planet. Sci. Lett., 277, 29-37. 

Leconte, J., Porget, F., Gharnay, B., Wordsworth, R., & Pottier, A. (2013). Increased insolation 
threshold for runaway greenhouse processes on Earth-like planets. Nature, 504, 268-271. 

Lenardic, A., Jellinek, M., & Moresi, L.-N. (2008). A climate change induced transition in the 
tectonic style of a terrestrial planet. Earth Planet. Sci. Lett., 271, 3442. 

Melosh, H. J., & Vickery, A. M. (1989). Impact erosion of the primordial atmosphere of Mars. 
Nature, 338, 487-489. 



- 52 - 


Menou, K. (2015). Climate stability of habitable Earth-like planets. Earth Planet. Sci. Lett., 4^9, 
20-24. 

Mills, B., Lenton, T. M., & Watson, A. J. (2014). Proterozoic oxygen rise linked to shifting balance 
between seafloor and terrestrial weathering. Proc. Natl. Aead. Sei., Ill, 9073-9078. 

Mills, B., Watson, A. J., Goldblatt, C., Boyle, R., &: Lenton, T. M. (2011). Timing of Neoproterozoic 
glaciations linked to transport-limited global weathering. Nature Geoseienee, 4, 861-864. 

Nakajima, S., Hayashi, Y.-Y., &: Abe, Y. (1992). A study on the ‘runaway greenhouse effect’ with 
a one-dimensional radiative-convective equilibrium model. J. Atmos. Sci., 49, 2256-2266. 

Newman, M. J., &: Rood, R. T. (1977). Implications of solar evolution for the earth’s early atmo¬ 
sphere. Science, 198, 1035-1037. 

Noack, L., & Breuer, D. (2014). Plate tectonics on rocky exoplanets: Influence of initial conditions 
and mantle rheology. Planet. Spaee Sci., 98, 41-49. 

Olson, P., Yuen, D. A., & Balsiger, D. (1984). Mixing of passive heterogeneities by mantle convec¬ 
tion. J. Geophys. Res., 89, 425-436. 

O’Neill, C., Jellinek, A. M., &: Lenardic, A. (2007). Conditions for the onset of plate tectonics on 
terrestrial planets and moons. Earth Planet. Sci. Lett., 261, 20-32. 

O’Neill, C., & Lenardic, A. (2007). Geological consequences of super-sized Earths. Geophys. Res. 
Lett, 34, 19204-19208. 

Pavlov, A. A., Kasting, J. F., Brown, L. L., Rages, K. A., & Freedman, R. (2000). Greenhouse 
warming by CH 4 in the atmosphere of early Earth. J. Geophys. Res., 105, 11981-11990. 

Riebe, C. S., Kirchner, J. W., & Finkel, R. C. (2004). Erosional and climatic effects on long-term 
chemical weathering rates in granitic landscapes spanning diverse climate regimes. Earth 
Planet. Sci. Lett., 224 , 547-562. 

Ronov, A. B., & Yaroshevsky, A. A. (1969). Chemical composition of the Earth’s crust. In 
P. Hart (Ed.), The Earth’s Grust and Upper Mantle (pp. 37-57). Washington DC: American 
Geophysical Union volume 13. 

Sagan, C., & Mullen, G. (1972). Earth and Mars: Evolution of Atmospheres and Surface Temper¬ 
atures. Science, 177, 52-56. 

Sleep, N. H., & Zahnle, K. (2001). Carbon dioxide cycling and implications for climate on ancient 
Earth. ,J. Geophys. Res., 106, 1373-1400. 


Sleep, N. H., Zahnle, K., &: Neuhoff, P. S. (2001). Initiation of clement surface conditions on the 
earliest Earth. Proe. Natl. Acad. Sci., 98, 3666-3672. 



- 53 - 


Stallard, R. F., & Edmond, J. M. (1983). Geochemistry of the Amazon: 2. The influence of geology 
and weathering environment on the dissolved load. J. Geophys. Res., 88, 9671-9688. 

Stamenkovic, V., Breuer, D., & Spohn, T. (2011). Thermal and transport properties of mantle rock 
at high pressure: Applications to super-earths. Icarus, 216, 572 - 596. 

Staudigel, H., Hart, S. R., Schmincke, H.-U., & Smith, B. M. (1989). Cretaceous ocean crust 
at DSDP Sites 417 and 418: Carbon uptake from weathering versus loss by magmatic 
outgassing. Geochim. Cosmochim. Acta, 53, 3091-3094. 

Stein, C., Lowman, J. P., &: Hansen, U. (2013). The influence of mantle internal heating on 
lithospheric mobility: Implications for super-Earths. Earth Planet. Sci. Lett., 361, 448-459. 

Strogatz, S. H. (1994). Nonlinear Dynamics and Chaos. Cambridge, MA: Westview Press. 

Tackley, P. (2000). The quest for self-consistent generation of plate tectonics in mantle convection 
models. In M. A. Richards, R. Gordon, &: R. van der Hilst (Eds.), History and Dynamics of 
Global Plate Motions, Geophys. Monogr. Ser. (pp. 47-72). Washington, DC: Am. Geophys. 
Union volume 121. 

Tackley, P. (2007). Mantle geochemical geodynamics. In D. Bercovici, & G. Schubert (chief edi¬ 
tor) (Eds.), Treatise on Geophysics (pp. 437-505). New York: Elsevier volume 7, Mantle 
Dynamics. 

Tajika, E. (1998). Climate change during the last 150 million years: reconstruction from a carbon 
cycle model. Earth Planet. Sci. Lett., 160, 695-707. 

Tajika, E., & Matsui, T. (1990). The evolution of the terrestrial environment. In H. E. Newsom, 
& J. H. Jones (Eds.), Origin of the Earth (pp. 347-370). 

Tajika, E., & Matsui, T. (1992). Evolution of terrestrial proto-C 02 atmosphere coupled with 
thermal history of the earth. Earth Planet. Sci. Lett., 113, 251-266. 

Takai, K., Nakamura, K., Toki, T., Tsunogai, U., Miyazaki, M., Miyazaki, J., Hirayama, H., Naka- 
gawa, S., Nunoura, T., & Horikoshi, K. (2008). Cell proliferation at 122°C and isotopically 
heavy CH 4 production by a hyperthermophilic methanogen under high-pressure cultivation. 
Proc. Natl. Acad. Sci., 105, 10949-10954. 

Trail, D., Watson, E. B., & Tailby, N. D. (2011). The oxidation state of Hadean magmas and 
implications for early Earth’s atmosphere. Nature, 4^0, 79-82. 

Turcotte, D., & Schubert, G. (2002). Geodynamics. (2nd ed.). New York: Cambridge Univ. Press. 

Valencia, D., & O’Connell, R. J. (2009). Convection scaling and subduction on Earth and super- 
Earths. Earth Planet. Sci. Letts., 286, 492-502. 



- 54 - 


Valencia, D., O’Connell, R. J., &: Sasselov, D. D. (2007). Inevitability of Plate Tectonics on Super- 
Earths. Astrophys. J., 670, L45-L48. 

van Heck, H. J., &: Tackley, P. J. (2011). Plate tectonics on super-Earths: Equally or more likely 
than on Earth. Earth Planet. Sci. Lett., 310, 252-261. 

Volk, T. (1987). Eeedbacks between weathering and atmospheric CO 2 over the last 100 million 
years. Am. J. Sei, 287, 763-779. 

Wade, J., & Wood, B. J. (2005). Core formation and the oxidation state of the Earth. Earth 
Planet. Sci. Lett., 236, 78-95. 

Walker, J., Hayes, P., & Hasting, J. (1981). A negative feedback mechanism for the long-term 
stabilization of Earth’s surface temperature. J. Geophys. Res., 86, 9776-9782. 

West, A. J. (2012). Thickness of the chemical weathering zone and implications for erosional and 
climatic drivers of weathering and for carbon-cycle feedbacks. Geology, 40 , 811-814. 

West, A. J., Galy, A., & Bickle, M. (2005). Tectonic and climatic controls on silicate weathering. 
Earth Planet. Sei. Lett., 235, 211-228. 

Willenbring, J. K., Codilean, A. T., &: McElroy, B. (2013). Earth is (mostly) flat: Apportionment 
of the flux of continental sediment over millennial time scales. Geology, 4 I, 343-346. 

Wolf, E. T., & Toon, O. B. (2013). Hospitable Archean Climates Simulated by a General Circulation 
Model. Astrobiology, 13, 656-673. 

Wolf, E. T., &: Toon, O. B. (2014). Delayed onset of runaway and moist greenhouse climates for 
Earth. Geophys. Res. Lett., 4 I, 167-172. 


This preprint was prepared with the AAS MI 5 X macros v5.2. 



